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from  the  field  campaign  of  the  Unified  Physical  Parameterization  for  Extended  Forecast 
(UPPEF  2012).  We  examined  the  general  characteristics  of  the  observed  stratocumulus- 
topped  boundary  layers  and  identified  cases  for  future  in-depth  studies.  We  first 
determined  the  boundary  layer  heights  from  all  sounding  profiles.  The  results  indicated 
the  sharp  westward  increase  of  boundary  layer  height  is  limited  to  ~200  km  offshore  with 
an  average  slope  of  2m  per  kilometer.  Substantial  west-east  spatial  variability  of 
thermodynamic  properties  is  also  observed  from  vertical  profiles.  The  sea  surface 
temperature  (SST)  in  the  region  varied  significantly.  Near  surface  measurements  over  the 
warm  and  cold  SST  regions  for  three  cases  were  analyzed  and  compared.  This  research 
found  that  the  small  scale  variability  in  the  SST  resulted  in  significant  variation  in  the 
surface  exchange  of  sensible  and  latent  heat  fluxes,  and  wind  stress.  Such  variability 
makes  it  difficult  to  correctly  parameterize  surface  fluxes.  The  presence  of  cool 
downdrafts  in  the  upper  cloud  layer  is  evident  in  the  joint  probability  density  distribution 
and  buoyancy  flux  profiles.  The  results  from  this  research  will  guide  future  cloud 
parameterization  scheme  development  in  forecast  models. 
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I.  INTRODUCTION 

A.  STRATOCUMULUS-TOPPED  BOUNDARY  LAYERS  AND  MILITARY 
OPERATIONS 

Stratocumulus  clouds,  derived  from  two  Latin  words  meaning  “layer”  and  “heap,” 
are  essentially  low,  individual  convective  elements  that  are  composed  in  a  layered  form 
under  a  shallow,  sharp  temperature  inversion.  Stratocumulus  topped  boundary  layers 
(STBL)  are  common  over  the  oceans  that  are  west  of  continents  on  the  eastern  part  of  a 
subtropical  high,  and  over  the  Arctic  Ocean  in  the  summer.  Stratocumulus  clouds  cover 
approximately  20%  of  the  Earth’s  surface  and  greatly  affect  energy  and  radiation 
balances  between  the  ocean  and  atmosphere.  They  are  strong  reflectors  of  incoming  solar 
radiation  and  weak  emitters  of  outgoing  longwave  radiation,  which  have  a  large  impact 
on  the  Earth’s  radiative  balance.  Recent  observations  show  that  low  clouds  reduce  the 
global  annual  radiation  balance  by  an  average  of  15  Wm'  (Klein  and  Hartmann  1993). 
There  are  few  areas  where  these  clouds  are  not  important  climato logically  (Wood  2012). 
The  ability  to  understand  the  complete  details  on  the  stratocumulus  cloud  processes  is 
very  difficult.  Much  research  has  developed  and  refined  the  theories  behind  the  STBL, 
including  an  original  proposal  by  Lilly  (1968)  known  as  the  mixed  layer  theory  (Brost  et 
al.  1982)  and  is  used  as  basis  for  other  scientists. 

Stratocumulus  clouds  have  been  studied  using  aircraft  for  the  past  ninety  years. 
The  first  aircraft  flight  that  researched  stratocumulus  clouds  took  place  on  September  10, 
1923  (Kloesel  1992).  This  flight  was  conducted  by  the  Naval  Air  Fleet  in  San  Diego,  CA 
with  the  goal  of  quantifying  and  forecasting  humidity  and  temperature  aloft.  The  data 
from  the  flight  became  a  starting  point  for  the  Navy  to  recognize  that  the  low 
stratocumulus  cloud  deck  which  primarily  existed  during  summer  months  had  great 
strategic  importance.  They  expended  many  man-hours  to  research  these  regions  of 
stratocumulus  clouds.  The  cloud  studies  generally  took  place  off  the  west  coast  of 
continents  with  cold  ocean  currents,  such  as  the  Californian  coast,  where  eighty  percent 
of  the  world’s  stratocumulus  clouds  exist  (Wood  2012).  Relatively  recent  extensive 
studies  on  stratocumulus  clouds  are  Fire  International  Satellite  Cloud  Climatology  Project 
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(ISCCP)  Regional  Experiment,  or  just  called  FIRE,  in  June- July  1987,  Atlantic 
Stratocumulus  Transition  Experiment  (ASTEX)  in  June  1992  (Kloesel  1992),  and  Second 
Dynamics  and  Chemistry  of  Marine  Stratocumulus  field  study  (DYCOMS-II)  in  July 
2001  (Stevens  et  al.  2003).  The  most  recent  experiment  is  Unified  Physical  Parameters 
for  Extended  Forecast  (UPPEF)  in  August- September  2012,  which  this  thesis  discusses  in 
detail.  The  past  research  has  given  a  solid  qualitative  understanding  on  formation, 
maintenance,  and  dissipation  of  stratocumulus  clouds.  However,  there  are  still  many 
uncertainties  and  in  order  to  accurately  parameterize  the  cloud,  additional  quantitative 
measurements  are  needed  to  model  the  net  effect  of  all  the  physical  processes  involved. 

In  recent  years,  various  research  groups  developed  PDF-related  cloud 
parameterizations  for  models,  mostly  using  larger  eddy  simulations.  However,  few  of  the 
parameterizations  are  evaluated  against  the  limited  aircraft  observations.  It  is  considered 
that  the  PDF  approach  may  lead  to  a  unified  cloud  parameterization  that  would  apply  to 
persistent  boundary  layer  clouds.  Also,  ensemble  forecasts  have  recently  become  state-of- 
the-art,  but  can  only  define  model  uncertainty  to  a  certain  extent.  It  is  known  that  most 
general  circulation  models  (GCM)  do  not  consistently  model  the  cloud  top  interface 
effectively  in  the  marine  stratocumulus  regions  (Stevens  2002).  Many  of  the  physical 
processes  in  stratocumulus  clouds  occur  on  time  and  spatial  scales  much  smaller  than  the 
GCM  time  step  and  grid  size.  Typically,  the  highest  resolution  numerical  weather 
prediction  (NWP)  models  have  horizontal  grid-box  scales  of  a  few  kilometers,  and  the 
spatial  scale  of  cloud  formation  and  dissipation  is  at  least  one  order  of  magnitude  smaller. 
This  shows  the  need  for  cloud  processes  to  be  parameterized  as  part  of  subgrid-scale 
processes.  Subgrid-scale  processes  are  defined  as  all  processes  that  cannot  be  resolved 
explicitly  with  present  or  future  models  (Kalnay  2003).  These  processes  include  turbulent 
mixing,  condensation,  evaporation,  friction,  and  radiation,  and  occur  on  scales  from  a  few 
meters  to  a  few  hundred  meters,  which  are  all  smaller  than  the  current  model  resolution. 
Ultimately,  enhanced  model  physics  and  realistic  lower  boundary  conditions,  such  as  sea 
surface  temperature,  will  improve  stratocumulus  cloud  and  other  subgrid-scale 
parameterizations . 
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Naval  forces  operate  in  many  regions,  whether  close  to  coast  in  littoral  regions  or 
over  the  open  ocean,  and  it  is  imperative  to  understand  the  battlespace  environment, 
especially  the  marine  boundary  layers.  The  warfighter  depends  on  accurate  models  to 
achieve  mission  success.  In  order  for  the  Navy  to  properly  deliver  a  complete  and 
representative  forecast  model  to  operational  forces,  the  boundary  layer  environment  must 
be  fully  known.  Most  naval  operations  are  conducted  in  the  boundary  layer — a  shallow 
but  dynamically  active  layer  of  the  atmosphere.  The  boundary  layer  is  also  the  layer  with 
persistent  turbulence  as  it  directly  interacts  with  the  Earth’s  surface  and  contributes  to 
most  daily  weather  phenomena,  such  as  precipitation,  fog,  significant  temperature 
changes  and  cloud  cover.  Many  naval  missions,  such  as  flight  operations,  sensor 
deployments  and  detection  of  electro-magnetic  signals,  communication,  and  ship 
tracking,  are  impacted  when  operating  under  solid  stratocumulus  clouds.  Remote  sensing, 
which  is  heavily  relied  on  for  many  tactical  military  operations,  is  directly  influenced  by 
the  presence  of  cloud  cover.  If  the  boundary  layer  characteristics  are  known,  relevant 
information  can  be  exploited  and  effectively  portrayed  to  ensure  mission  success. 

B.  THESIS  OBJECTIVES 

During  September  2012,  extensive  aircraft  measurements  were  made  off  the  coast 
of  Monterey,  CA  as  part  of  UPPEF  2012,  which  the  Office  of  Naval  Research  (ONR) 
sponsored  as  a  Departmental  Research  Initiative  (DRI).  UPPEF  2012  field  measurement 
utilized  the  Twin  Otter  research  aircraft  operated  by  the  Naval  Postgraduate  School 
(NPS)  Center  for  Interdisciplinary  Remotely-Piloted  Aircraft  Studies  (CIRPAS)  for 
twelve  research  flights,  ten  of  which  were  made  in  cloudy  boundary  layers  with  various 
cloud  fractions.  In  addition  to  NPS,  the  Naval  Research  Laboratory  (NRL)  in  Monterey, 
CA  and  University  of  California,  Irvine  (UCI)  also  participated  in  UPPEF  2012.  The  goal 
of  the  experiment  was  to  study  the  cloud  structure,  turbulence,  and  other  related  physical 
processes  in  broken  stratocumulus  using  aircraft-based  instruments  and  platforms. 

This  thesis  work  analyzes  the  UPPEF  2012  dataset  in  order  to  seek  answers  to 
some  of  the  scientific  questions  listed  in  the  earlier  section.  The  ultimate  goal  is  to  help 
improve  cloud  parameterizations  in  forecast  models.  This  research  provides  the  initial 

analyses  of  the  aircraft  measurements  from  UPPEF  2012.  The  first  step  is  to  examine  all 
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measurements  and  identify  specific  cases  for  focused  in-depth  analyses.  A  variety  of  data 
analyses  methods  is  used  including  statistical  analyses  and  spectral  analyses,  along  with 
different  graphical  presentations  are  used  to  depict  the  vertical  and  horizontal  variability 
of  the  cloud-topped  marine  boundary  layer. 

This  thesis  is  structured  as  following.  It  begins  with  an  introduction  on  the 
atmospheric  boundary  layer  and  physical  processes  within  stratocumulus  clouds.  An 
overview  of  UPPEF  2012  is  presented  next,  followed  by  an  analysis  of  the  data  collected 
and  used  for  this  thesis  research.  Then,  two  case  studies  are  presented  on  boundary  layer 
variability  and  three  cases  are  discussed  showing  the  sea  surface  temperature  variability. 
Next,  the  turbulence  characteristics  are  analyzed  in  cloud  and  at  low  levels.  Finally, 
conclusions  are  made  with  recommendations  for  future  studies. 
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II.  BACKGROUND 


A.  THE  ATMOSPHERIC  BOUNDARY  LAYER 

The  lowest  portion  of  the  atmosphere  is  considered  the  atmospheric  boundary 
layer  (ABL),  also  known  as  the  planetary  boundary  layer  (PBL)  or  simply  the  boundary 
layer  (BL).  The  BL  is  “part  of  the  troposphere  that  is  directly  influenced  by  the  presence 
of  the  earth’s  surface,  and  responds  to  surface  forcing  with  a  timescale  of  about  an  hour 
or  less”  (Stull  1988).  It  is  extremely  important  to  the  climate,  weather  and  air  quality. 
However,  the  structure  of  the  BL  is  often  very  complex.  The  BL  height  characterizes  the 
vertical  extent  of  mixing  in  the  boundary  layer  and  the  level  where  momentum,  water, 
and  trace  substances  are  exchanged  with  the  free  troposphere.  The  bottom  boundary  of 
the  atmosphere  is  the  Earth’s  surface.  Atop  the  BL  is  a  capping  inversion,  a  stable  layer 
between  the  BL  below  and  the  free  atmosphere  above,  which  prevents  turbulence  mixing 
into  the  deeper  troposphere.  Exchanges  of  pollutants  and  moisture  with  the  free 
troposphere  are  conducted  through  entrainment  at  the  BL  top  (Hobbs  and  Wallace  2006). 
Over  the  ocean,  the  BL  is  often  referred  to  as  the  marine  boundary  layer  (MBL).  MBL 
tends  to  be  cloud-topped  with  stratocumulus,  cumulus,  and  stratus  clouds  (Hennemuth 
and  Lammert  2005). 

The  BL  height  is  a  key  parameter  used  in  many  models.  There  are  many  methods 
used  to  determine  the  BL  height.  Rawinsondes  are  commonly  used  to  determine  the  BL 
height,  since  routine  soundings  are  done  worldwide  with  high  accuracy  and  relatively  low 
cost.  Rawinsondes  observe  only  on  a  single  line  (Hennemuth  and  Lammert  2005)  but 
offer  vertical  profiles  of  temperature  and  humidity  through  the  inversion  and  up  to  tens  of 
kilometers  in  regions  mostly  inaccessible  to  humans  (Flores  et  al.  2013).  The  BL  height 
estimation  for  UPPEF  2012  uses  primarily  rawinsondes  and  aircraft  measurements. 
Aircraft  provide  concurrent  measurements  of  many  variables  that  can  be  used  to  estimate 
the  height  (Seibert  2000).  They  are  also  well-suited  to  measure  mean  and  turbulence 
quantities  along  horizontal  legs  at  various  height  levels  (Brost  et  al.  1982).  The  BL  height 
is  constantly  changing  both  spatially  and  temporally  depending  on  the  geographic 
location,  topography,  surface  cover,  time  of  day,  season,  and  weather  patterns.  The  BL 
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can  vertically  extend  from  tens  of  meters  to  a  few  kilometers  in  various  regions  of  the 
world  (Stull  1988).  It  can  even  collapse  and  turn  “upside  down”  in  extremely  stable 
conditions  (Sun  et  al.  2002,  2004).  In  case  of  the  STBL,  the  BL  height  is  consistent  with 
the  cloud  top  height  below  the  strong  capping  inversion  (Hennemuth  and  Lammert  2005). 

B.  CONCEPTUAL  MODEL  OF  THE  CALIFORNIA  MARINE  BOUNDARY 

LAYER 

The  MBL  varies  significantly  in  the  coastal  region  due  to  the  complex  interaction 
among  the  coastal  topography,  air-sea  temperature  and  surface  roughness  contrast,  and 
synoptic  conditions.  California  is  topographically  diverse  with  mountains  and  valleys.  To 
the  west  of  the  Sierra  Nevada  mountain  range  and  the  Central  Valley  lie  the  Coastal 
Ranges  along  the  northern  and  central  CA  coastline.  These  include  the  Diablo  Range,  just 
east  of  San  Francisco  and  the  Santa  Cruz  Mountains  located  to  the  south  of  San 
Francisco.  These  mountain  ranges  are  home  to  the  giant  redwoods,  the  tallest  trees  on  the 
planet.  The  coastal  mountain  ranges  affect  the  local  surface  wind  flow  and  the  boundary 
layer  characteristics  off  the  CA  coast. 

Beardsley  et  al.  (1987)  proposed  a  conceptual  model  meant  to  describe  the  mean 
summertime  state  of  the  lower  atmosphere  over  the  Eastern  Pacific  (EPAC)  during 
periods  of  persistent  northwesterly  winds,  shown  in  Figure  1.  The  model  assumes  a 
straight  coast,  with  an  adjacent  closed  coastal  mountain  range,  to  account  for  the  Coastal 
Range.  The  domain  is  divided  into  three  regions  parallel  to  the  coast:  a  nearshore  zone 
(0-20km),  an  intermediate  zone  (20-200km),  and  an  offshore  zone  (200-1 500km).  In  the 
offshore  zone,  the  marine  inversion  tilts  slowly  down  toward  the  coast.  Furthest  west, 
near  1500km  offshore,  the  inversion  is  at  its  maximum  at  1500m  and  decreases  to  300m 
as  the  intermediate  zone  is  approached.  The  eastward  tilt  reflects  the  effect  of  the 
subsidence  from  the  summertime  EPAC  high  and  north  to  northwesterly  surface  winds 
with  speeds  between  5-7.5  ms'1.  Stratus-type  clouds  and  fog  only  occur  below  the 
inversion.  In  the  intermediate  zone,  the  surface  wind  is  typically  faster,  about  10  ms'1,  off 
the  northern  CA  coast  than  anywhere  else  along  the  U.  S.  west  coast.  The  maximum  wind 
speed  in  this  area  is  partly  due  to  the  EPAC  high  located  directly  to  the  west  of  northern 
CA.  Other  factors  are  towards  the  Gulf  of  Alaska,  synoptic  features  traveling  to  the  east 
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often  interrupt  the  tendency  for  the  surface  isobars  to  parallel  the  coast,  and  towards  the 
south,  as  the  CA  coast  turns  eastward,  large-scale  divergence  increases  within  the  marine 
layer.  The  effect  contributes  to  the  maximum  alongshore  flow  along  the  northern  CA 
coast.  There  is  also  a  diurnal  variation  in  the  wind,  by  1  ms'1,  with  the  strongest  northerly 
winds  in  the  morning.  Higher  wind  speeds  in  the  intermediate  zone  imply  that  the  marine 
layer  must  tilt  more  steeply  toward  the  east  by  geostrophy. 

1000 m 
A.  200 


Figure  1.  Conceptual  model  of  summertime  MBL  over  northeastern  PAC  during 
periods  of  persistent  northerly  winds  (From  Beardsley  et  al.  1987). 

The  nearshore  zone  is  characterized  by  land-sea  diurnal  variability  (Beardsley  et 

al.  1987).  A  separate  conceptual  model  is  shown  in  Figure  2.  The  maximum  variability 

occurs  over  the  beach  and  decreases  rapidly  offshore.  Over  the  narrow  coastal  plain, 

nighttime  infrared  (IR)  cools  and  stabilizes  the  lower  atmosphere,  resulting  in  calm  winds 

in  the  lower  50m  over  land.  At  night,  the  northwesterly  winds  weakly  increase  with 

height,  from  50  m  to  the  inversion,  without  any  significant  cross-coast  flow.  In  the  early 

morning,  there  is  a  sharp  and  level  inversion  and  the  maximum  northwesterly  wind  is 

20km  offshore.  After  sunrise,  the  sun  heats  the  surface  and  destabilizes  the  coastal  air 

causing  a  weak  density-driven  cross-coast  circulation  to  set  in  and  the  maximum 

alongshore  wind  to  penetrate  onshore.  There  is  a  weak  return  flow  above  the  eroded 
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inversion  and  the  subsidence  branch  of  this  circulation  cause  the  inversion  to  decrease  to 
a  minimum  height  of  50-100  m  over  the  beach  (Beardsley  et  al.  1987).  Inland,  turbulent 
mixing  and  surface  heating  causes  the  inversion  to  lift  slowly  and  diffuse.  The  maximum 
wind  speed  is  located  just  below  the  base  of  the  inversion  offshore,  while  it  is  broadly 
located  between  the  inversion  base  and  the  ground  inland.  As  the  sun  goes  down,  the 
stability  over  land  increases  due  to  radiative  cooling  and  the  weak  cross  flow  ceases. 
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Figure  2.  Conceptual  model  of  variability  of  MBL  height  and  low-level  wind  in  the 
nearshore  zone  during  (a)  night  and  (b)  day  (From  Beardsley  et  al.  1987). 

C.  SEPTEMBER  CLIMATOLOGY  OFF  NORTHERN  CALIFORNIA 
1.  Synoptic  Pattern 

The  northeastern  Pacific  Ocean  is  dominated  by  a  strong  Eastern  Pacific  (EPAC) 
high  during  summer  to  early  fall  interacting  with  a  heat  low  over  the  deserts  of  the 
southwestern  United  States  causing  a  strong  and  persistent  pressure  gradient  at  the 
surface.  The  pressure  gradient  drives  a  steady  northwesterly  surface  wind  along  the 
California  coast.  Figure  3a  illustrates  the  climatological  sea  level  pressure  (SLP)  and 
surface  synoptic  patterns:  EPAC  high,  thermal  low  over  the  desert  states,  and  the 
resulting  surface  pressure  gradient.  The  surface  wind  for  September  is  shown  in  Figure 
3b,  illustrating  the  strongest  northwesterly  wind  at  about  400  km  offshore.  These  figures 
were  generated  using  the  National  Weather  Service  (NWS)  National  Centers  for 
Environmental  Prediction  and  the  National  Science  Foundation  (NSF)  National  Center 
for  Atmospheric  Research  (NCEP/NCAR)  Reanalysis  database  located  in  the  National 
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Oceanographic  &  Atmospheric  Administration  (NOAA)  Physical  Science  Division  and 
represents  the  average  conditions  for  September  with  almost  30  year  of  data,  from  198 1— 
2010.  Brost  et  al.  (1982)  notes  this  pressure  gradient  is  probably  the  strongest  and  most 
persistent  summertime  pressure  gradient  observed  in  the  entire  North  America. 


NCEP/NGAR  RaonalyBia  NCEP/NCAR  RaonolyBia 

San  Level  Pressure  (mb)  Climutology  1981-2010  cllmo  Surfoce  Vector  Wind  (m/s)  Climatology  19B1-2010  climo 


Figure  3.  (a)  Sea  level  pressure  and  (b)  surface  wind  speed  and  wind  vector 

climatology  for  September,  generated  from  NCEP/NCAR  reanalysis  dataset. 

The  mean  summer  atmospheric  regime  is  generally  interrupted  by  five  anomalous 
synoptic  or  mesoscale  events  (Beardsley  et  al.  1987).  First,  a  strengthened  gradient  with 
enhanced  northwesterly  flow  can  occur  for  several  days.  This  causes  a  stronger  than 
normal,  upwelling  favorable  wind  event.  Second,  a  weak,  quick  moving  upper  level 
trough  passes  and  produces  stronger  than  normal  northwesterly,  upwelling  favorable 
winds,  but  limited  to  a  day  or  so.  Third,  an  internal  Kelvin  wave  will  produce  a  steady 
poleward  progression  of  southerly  winds,  which  is  considered  a  wind  relaxation  or  wind 
reversal  event.  The  Kelvin  wave  causes  the  marine  layer  to  rise  along  the  coast  and  out  to 
50  km  offshore.  Satellite  pictures  show  overcast  stratus  to  mark  the  leading  edge  of  the 
lifted  marine  layer  with  a  clear  area  to  the  north.  The  clear  area  contains  strong  northerly 
winds  and  a  lower  marine  layer.  Kelvin  waves  typically  travel  6  ms'1.  As  the  leading  edge 
passes,  the  surface  winds  suddenly  switch  to  southerly  allowing  overcast  stratus  to  set  in 
along  the  coast.  A  few  days  later,  the  Kelvin  wave  passes  and  winds  suddenly  reverse 
back  to  northerly.  Fourth,  another  wind  reversal  event  is  caused  by  alongshelf  southerly 
winds  developing  inshore  along  the  northern  CA  coast.  This  event  is  contributed  to 
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coastally  trapped  gravity  currents  that  form  at  the  southern  CA  coast  and  surge 
northward.  This  wind  reversal  is  limited  to  a  coastal  zone  of  less  than  100km  wide.  Along 
the  coast,  the  inversion  rises  in  spurts  with  the  onset  of  the  southerly  winds  and  overcast 
stratus  develops  in  the  raised  inversion.  When  the  event  ends,  the  inversion  height 
decreases  and  the  wind  returns  to  normal.  Finally,  the  smallest  scale  wind  reversal  occurs 
for  a  few  hours  in  the  morning  and  only  in  the  very  nearshore  zone. 

2.  Sea  Surface  Temperature  (SST)  Characteristics 

SST  strongly  affects  the  earth’s  climate  and  radiation  budget,  and  influences 
cloud  cover  and  turbulent  heat  fluxes.  California  is  known  to  have  relatively  cold  coastal 
waters.  The  combined  Coriolis  effect  and  the  wind  stress  from  the  maximum 
northwesterly  alongshore  wind  off  central  and  northern  CA  pushes  the  upper  layer  of 
water  offshore  to  the  southwest.  The  upper  layer  of  water  is  replaced  by  cold  bottom 
water  through  upwelling,  resulting  in  low  sea  surface  temperatures  immediately 
offshore.  The  marine  atmospheric  boundary  layer  (MABL)  over  the  cool  SST  region  is 
generally  colder  and  prone  to  saturation,  especially  between  Point  Conception  and  Cape 
Blanco  (Beardsley  et  al.  1987).  The  combination  of  increased  subsidence  from  the  strong 
EPAC  high  and  a  cool  MABL  creates  a  very  strong  inversion  at  the  MABL  top  during  the 
summertime,  which  in  turn  traps  the  cool,  moist  air  close  to  the  surface,  causing  much  of 
the  Californian  coastline  to  be  covered  with  thick  stratus  or  fog  for  a  majority  of  the  day. 

Climatologically,  during  September,  SST  is  lowest  just  off  the  CA  coast  and 
increases  to  the  southwest,  as  depicted  in  Figure  4,  generated  using  NCEP/NCAR 
reanalysis  database.  The  upwelled  water  near  the  coast  is  much  colder  (less  than  15  °C) 
than  the  water  offshore  and  sharp  ocean  fronts  often  occur  as  a  boundary  between  the 
warm  and  cold  water  masses  as  seen  from  daily  SST  fields. 
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NQAA  Extends d  5 ST 


Figure  4.  September  SST  climatology  using  NCEP/NCAR  reanalysis  database. 

Filaments,  also  known  as  squirts,  plumes,  and  tongues,  can  also  occur  and  affect 
the  BL.  Upwelling  filaments  originate  on  the  CA  shelf  and  are  narrow  jets  that  advect 
cold  coastal  water  several  hundreds  of  kilometers  offshore  (Flament  et  al.  1985).  There 
are  three  types  of  water  across  a  filament:  high-salinity  warm  water  to  the  south,  low- 
salinity  warmer  water  to  the  north,  and  cold  lower-salinity  water  in  the  filament’s  axis 
(Flament  et  al.  1985).  The  filament’s  boundaries  are  asymmetric.  The  southern  boundary 
of  cold  water  is  constrained  by  a  sharp  front  less  than  1km  wide  and  the  northern 
boundary  has  a  gradual  transition  from  cold  to  warm  water  over  20km.  Instabilities  and 
strong  turbulence  can  occur  over  the  sharp  front. 

There  are  extensive  sheets  of  stratocumulus  in  relatively  shallow,  well-mixed  BL 
over  cold  SST.  SST  is  one  of  the  most  important  variables  influencing  cloud  fraction  in 
marine  stratocumulus.  It  influences  the  stratocumulus  clouds  in  two  time  scales  (Hegg  et 
al.  2004).  On  short  time  scales,  higher  SSTs  are  associated  with  higher  heat  and  water 
vapor  fluxes  that  moisten  the  MBL.  They  also  produce  warmer  temperatures  and  greater 
water  vapor  mixing  ratios  at  the  cloud  base,  resulting  in  more  liquid  water  content  (LWC) 
and  thicker,  moister  clouds.  BL  cloudiness  is  correlated  with  upstream  SST  on  a 
timescale  of  approximately  a  day  (Pincus  et  al.  1997).  On  longer  time  scales,  higher  SST 
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initiates  the  stratocumulus  to  cumulus  transition  and  breaks  up  the  cloud  deck.  Diurnal 
variations  introduced  by  solar  radiation  can  decouple  the  cloud  layer  from  the  surface 
during  the  day  in  these  regions.  The  gradual  increase  of  SST  on  a  large  scale,  e.g.,  from 
the  CA  coast  southwest  toward  Hawaii,  deepens  the  boundary  layer  and  eventually  leads 
to  permanent  decoupling  of  the  cloud  layer  from  the  surface-based  mixed  layer  (SML) 
(e.g.,  Pincus  et  al.  1997,  Bretherton  et  al.  1995).  Moisture  builds  up  in  the  SML  which 
supports  cumulus  cloud  growth  from  the  top  of  the  SML.  Cumulus  clouds  either  help 
maintain  stratocumulus  by  supplying  it  with  moisture  from  the  SML  or  dissipate  by 
enhancing  the  entrainment  into  the  cloud  layer  (Wang  and  Lenschow  1995,  Bretherton  et 
al.  1995).  A  schematic  diagram  of  the  transition  of  stratocumulus  to  cumulus  cloud  is 
provided  in  Figure  5.  UPPEF  2012  was  primarily  focused  on  the  stratocumulus  cloud  in 
the  shallow,  well-mixed  layer  over  cold  SSTs  off  the  CA  coast. 


Figure  5.  An  illustration  of  the  transition  from  stratocumulus  to  trade  wind  cumulus 

(From  Bretherton  et  al.  1995). 


3.  Inversion  Characteristics 

There  is  a  strong  temperature  inversion  over  most  of  the  EPAC  that  clearly  marks 
the  boundary  between  subsiding  dry,  warm  air  above  and  cool,  moist  air  in  the  turbulent, 
marine  layer  below.  The  inversion  is  in  general  thickest  and  has  the  greatest  temperature 
difference  between  the  top  and  base  just  off  northern  CA  (Beardsley  et  al.  1987). 
Typically,  the  inversion  is  about  400  m  over  North  and  central  CA  coast,  rising  quickly  to 
the  west  to  a  maximum  of  about  2000  m  near  Hawaii  (Beardsley  et  al.  1987).  According 
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to  Brost  et  al.  (1982),  near  the  CA  shore  the  BL  is  shallow,  stably  stratified,  and  foggy. 
Father  out,  it  becomes  deeper,  near-neutral,  and  stratus  or  stratocumulus  topped.  Farther 
still,  toward  Hawaii  (HI),  the  cloud  layer  breaks  up  into  scattered  cumulus  (Figure  5). 
There  is  a  north-south  change  in  inversion  height  due  largely  to  entrainment  and 
mesoscale  vertical  motion.  The  east-west  change  in  inversion  slope  is  due  to  the 
average  sea  surface  temperature  increase  to  the  west  and  weaker  subsidence  away  from 
the  center  of  the  subtropical  high.  The  inversion  slope  can  cause  large  vertical  wind 
shears  (5-10  ms'1)  with  winds  being  faster  below  the  inversion  than  immediately  above 
(Beardsley  et  al.  1987).  The  relatively  low  inversion  in  the  stratocumulus-topped  MABL 
forces  the  wind  direction  to  be  constant  with  height  in  the  BL,  with  a  jump  in  wind 
direction  across  the  inversion  (Brost  et  al.  1982).  At  the  inversion,  wind  shear  is  large  and 
produces  turbulence  and  entrainment. 

D.  PHYSICAL  PROCESSES  IN  THE  STRATOCUMULUS-TOPPED 

BOUNDARY  LAYER 

1.  General  Characteristics 

Low  stratocumulus  clouds  are  found  primarily  over  the  oceans  and  on  the  eastern 
side  of  the  semi-permanent  subtropical  high  pressure  systems.  Stratocumulus  clouds  form 
just  below  a  strong  temperature  inversion  below  1500  m  in  altitude  (Wood  2012),  which 
marks  the  marine  boundary  layer  top  over  large  areas  of  the  oceans.  They  are  also 
normally  referred  to  as  low-level  clouds  primarily  driven  by  convective  instability  due  to 
radiative  cooling  at  the  cloud  top.  The  maximum  cloud  amount  occurs  during  the  summer 
months  (June,  July,  and  August)  when  the  low-level  thermal  stability  also  peaks  (Klein 
and  Hartmann  1993).  The  clouds  dissipate  by  processes,  such  as  entrainment  drying, 
decoupling,  or  drizzle  or  by  transitioning  to  another  cloud  type  (cumulus  for  example) 
(Wood  2012). 

The  vertical  structure  of  the  stratocumulus-topped  boundary  layer  (STBL)  is 
strongly  tied  to  the  associated  stratocumulus  clouds.  Fractional  cloud  coverage  is  greatest 
when  the  STBL  is  moderately  shallow  (height  is  500-1 000m)  and  well  mixed.  Figure  6 
shows  typical  vertical  thermodynamic  structures  of  water  vapor  ( q )  and  liquid  water 
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mixing  ratios  ( qi ),  equivalent  potential  temperature  (6e)  and  temperature  (7)  for 
summertime  shallow,  well-mixed  STBL  corresponding  to  the  photograph  of  solid 
stratocumulus  clouds  in  Figure  7.  The  cloud  top  corresponds  to  the  level  that  qi  decreases 
abruptly  from  its  maximum  with  a  strong  T  and  6e  inversion,  and  rapid  decrease  in  q.  The 
liquid  water  mixing  ratio  typically  increases  with  height  at  a  quasi-linear  rate  in 
stratocumulus  layers,  and  is  primarily  affected  by  temperature  and  weakly  dependent  on 
pressure  (Wood  2012).  The  boundary  layer  is  below  the  shallow  and  sharp  inversion 
normally  less  than  100m  thick. 


Figure  6.  Typical  thermodynamic  vertical  profiles  of  q,  qi,  0e,  and  T  for  summertime 

shallow,  well-mixed  STBL.  Dots  depict  horizontal  leg  means  from  constant 
level  measurements.  Significant  layers  and  boundaries  are  labeled 
(After  Nicholls  1984). 


The  STBL  is  capped  by  a  shallow  layer  called  the  entrainment  interfacial  layer 
(EIL)  which  starts  at  the  base  of  the  inversion  layer  atop  the  cloud.  This  layer  has  strong 
gradients  in  temperature,  humidity,  and  cloud  liquid  water.  The  universal  defining  feature 
of  the  EIL  is  the  strong  temperature  increase  (inversion)  and  moisture  decrease  with 
height  (Figure  6).  Virtual  potential  temperature  ( 6V )  and  total  water  content  (qr)  maybe 
strongly  correlated  in  the  EIL  (Nicholls  and  Leighton  1986).  In  Nicholls  and  Leighton 
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(1986),  the  EIL  is  often  subsaturated,  there  is  a  net  negative  buoyancy  flux,  and  vertical 
velocity  fluctuations  are  relatively  small  in  this  layer.  However,  the  entrainment  induced 
evaporative  cooling  may  also  result  in  strong  cool  downdrafts  as  illustrated  in  Wang  and 
Albrecht  (1994).  The  EIL  is  on  the  order  of  a  few  meters  to  tens  of  meters  thick,  but  is 
highly  variable  within  the  same  cloud.  The  difference  between  air  properties  at  the  top  of 
the  cloud  layer  and  the  EIL  are  usually  so  marked  and  the  distance  between  the  two 
layers  is  small,  it  is  often  represented  by  a  discontinuity  or  jump.  Nicholls  and  Leighton 
(1986)  noticed  that  data  from  flight  legs  close  to  the  cloud  top  had  quite  different 
characteristics  than  those  found  only  a  few  tens  of  meters  deeper  in  the  cloud.  Beneath 
the  EIL,  at  a  depth  of  about  50  m,  negatively  buoyant  downdrafts  are  established,  9V  and 
qr  are  often  positively  correlated,  and  vertical  velocity  variance  increases.  This  suggests 
radiatively  cooled  downdrafts  have  combined  some  drier  air  from  aloft  during  cloud  top 
entrainment. 


Figure  7.  Photograph  of  a  solid  stratocumulus  cloud  top  (From  Wood  2012). 

Low  level  marine  stratocumulus  clouds  are  important  contributors  to  the  earth’s 

radiation  budget.  Stratocumulus  clouds  have  a  high  albedo  (strongly  reflective)  compared 

to  the  ocean.  Areas  with  persistent  stratocumulus,  such  as  the  California  coast,  absorb 

twenty  percent  less  solar  energy  annually  than  other  places  at  the  same  latitude  (Hanson 

and  Gruber  1982).  Liquid  water  content  (LWC)  determines  the  cloud’s  optical  properties 
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and  is  a  critical  link  between  cloud  dynamics  and  radiative  effects.  As  optical  depth 
increases,  more  solar  radiation  is  reflected  to  space,  cooling  the  atmosphere  and  the 
surface.  Pincus  et  al.  (1997)  showed  that  clouds  with  large  optical  thickness  values  in  the 
morning  time  are  unlikely  to  break  up  during  the  afternoon. 

On  any  given  day,  there  can  be  large  and  mesoscale  variations  in  the 
stratocumulus  cloud  clearly  visible  in  satellite  imagery.  One  variation  from  solid  cloud 
off  the  CA  coast  is  a  rift  area.  Rifts  are  areas  with  low  reflectivity,  distinctly  cellular 
conditions  seen  from  satellite,  and  advect  with  the  BL  winds  (Sharon  et  al.  2006).  Rift 
regions  contain  broken  stratocumulus  clouds  surrounded  by  solid  stratocumulus  clouds. 
Another  characteristic  is  that  ship  tracks  have  enhanced  reflectivity  in  rift  areas  due  to  the 
increased  aerosols  from  the  ship  exhaust  and  a  decrease  in  drizzle.  This  enhanced 
reflectivity  of  ship  tracks  makes  it  easier  to  depict  rifts  in  satellite  images. 

2.  Key  Physical  Processes 

The  key  physical  processes  that  occur  in  the  STBL  are  longwave  cooling,  solar 
heating,  entrainment,  turbulent  mixing,  latent  heating,  evaporative  cooling,  surface  fluxes 
and  wind  shear,  illustrated  in  Figure  8.  Stratocumulus,  a  shallow  convective  cloud 
system,  is  maintained  by  the  generation  of  buoyant  instability  by  radiative  cooling  at  the 
cloud  top,  which  drives  turbulent  mixing  to  result  in  the  observed  temperature  and 
humidity  profiles,  causes  cloud  top  entrainment,  and  affects  surface  fluxes. 


Longwave 
cooling  Solar 


Figure  8.  Schematic  showing  the  key  processes  that  occur  in  the  STBL 

(From  Wood  2012). 
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a. 


Longwave  Cooling 


Longwave  cooling  is  the  leading  term  in  the  STBL  energy  budget  (Wood 
2012),  provides  a  major  source  of  buoyancy  that  drives  convective  instability  to  generate 
and  maintain  turbulence  in  the  cloud  layer  and  below  (Moeng  2000,  Wood  2012).  Most 
stratocumulus  clouds  contain  a  sufficiently  abundant  amount  of  liquid  water  so  they  are 
essentially  opaque  to  longwave  radiation.  The  longwave  volume  absorption  coefficient 
increases  quasi-linearly  with  qi.  Normally,  the  depth  that  the  longwave  radiation 
penetrates  in  stratocumulus  cloud  tops,  with  qj> 0.2  g  kg'1,  is  only  a  few  meters  to  tens  of 
meters.  Stratocumulus  clouds  tend  to  occur  under  a  dry  free  troposphere  with 
downwelling  longwave  radiative  fluxes  just  above  the  cloud  top  is  less  than  the  upwelling 
flux,  by  several  tens  to  a  little  over  100  W  m~2  (Wood  2012).  No  more  than  a  few  tens  of 
meters  down  from  the  cloud  top,  the  downwelling  and  upwelling  fluxes  are 
approximately  equal.  Thus,  the  upper  few  meters  of  the  cloud  experience  strong 
longwave  cooling.  The  average  flux  difference  in  this  layer  is  50-90  W  m'2  (Nicholls  and 
Leighton  1986,  Wood  2012).  Also,  near  the  cloud  base,  net  longwave  flux  convergence 
causes  weak  heating.  Away  from  the  cloud  layer,  the  longwave  flux  divergence  is  small 
and  does  not  significantly  contribute  to  the  net  heating  or  cooling  of  the  STBL.  In 
addition,  longwave  cooling  is  considered  a  key  dynamical  element  that  drives  overturning 
of  convective  circulations  in  the  clouds. 

b.  Solar  Heating 

Solar  heating,  also  called  solar  absorption,  is  another  major  component  in 
the  STBL  energy  budget  (Wood  2012)  and  is  the  primary  driver  of  the  STBL  diurnal 
cycle  (Turton  and  Nicholls  1987,  Wood  2012).  The  majority  of  solar  absorption  in  cloud 
occurs  in  the  near-infrared  band  of  the  electromagnetic  spectrum.  The  fraction  of  the 
absorbed  incoming  solar  radiation  by  stratocumulus  can  reach  up  to  15%  (Wood  2012). 
The  fractional  absorption  decreases  with  the  solar  zenith  angle  and  increases  with  cloud 
optical  depth  (to  a  limit).  Water  droplet  absorbs  close  to  half  of  the  solar  radiation  in  the 
cloud  layer.  For  all  cloud  optical  depths  and  solar  zenith  angles,  solar  heating  is  largest  at 
the  cloud  top  due  to  strong  scattering  of  solar  radiation  at  the  cloud  top  limits  absorption 
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lower  down.  In  most  stratocumulus  clouds  that  the  liquid  water  content  increases  with 
height,  solar  heating  decreases  much  slower  down  through  the  cloud  layer  than  longwave 
cooling.  Even  in  cases  of  strong  solar  insolation,  the  net  effect  of  radiation  acts  to 
destabilize  the  cloud  layer  (Wood  2012)  causing  continued  entrainment  as  the  layer 
warms  and  decouples  from  the  surface  (Pincus  et  al.  1997). 

c.  Turbulent  Mixing 

Turbulent  mixing  homogenizes  the  boundary  layer  and  tends  to  couple  the 
layer  to  the  surface  which  maintains  the  moisture  supply  to  the  cloud.  The  STBL’s  mean 
state  is  predominately  shaped  by  turbulent  fluxes  of  energy,  water  vapor,  and  liquid 
water.  The  turbulent  kinetic  energy  (TKE)  influences  the  entrainment  rate  of  the  free- 
tropospheric  air  into  the  STBL.  Profiles  of  vertical  turbulent  fluxes  are  important  to 
determine  STBL  properties.  A  classic  depiction  of  the  STBL  is  described  by  Lilly  in 
1968  as  the  mixed  layer  theory.  Lilly’s  theory  describes  the  vertical  structure  of  fluxes 
necessary  for  maintaining  a  well-mixed  layer.  For  a  completely  well-mixed  layer,  the 
vertical  energy  and  moisture  fluxes  are  linear  functions  of  height.  For  non-precipitation 
mixed  layers,  the  vertical  turbulent  flux  of  total  water  is  also  linear  with  height  (Wood 
2012).  The  primary  generator  of  TKE  in  the  STBL  is  the  buoyancy  flux,  which  is  a 
maximum  in  the  cloud  layer.  TKE  is  dissipative  and  not  conserved.  The  boundary  layer 
can  be  turbulent  only  if  there  are  specific  physical  processes  generating  the  turbulence. 
The  TKE  energy  balance  shows  that  buoyancy  production  is  largest  just  below  the  cloud 
top  and  near  the  surface  (Nicholls  1984).  TKE  is  transported  downward  from  the  cloud 
top  through  turbulent  mixing. 

Buoyancy  fluxes  are  caused  by  radiative  cooling  and  enhanced  by  latent 
heating.  In  the  STBL,  latent  heat  release  causes  the  buoyancy  flux  to  increase  just  above 
the  cloud  base.  The  maximum  buoyancy  flux  is  located  at  the  upper  part  of  the  cloud 
layer,  implying  that  cloud  top  cooling  is  a  dominant  mechanism  in  buoyancy  production. 
In  lower  levels  of  the  BL,  the  buoyancy  flux  decreases  and  is  small  and  can  be  slightly 
negative  (Nicholls  1989).  Radiative  cooling  at  the  cloud  top  results  in  negatively  buoyant 
eddies  and  entrains  drier  air  from  above  the  boundary  layer  as  they  become  downdrafts. 
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Negatively  buoyant  downdrafts  are  the  primary  convective  elements  in  a  cloud  layer. 
Updrafts  are  warmer  and  have  higher  LWC  than  the  downdrafts. 

d.  Entrainment 

Entrainment  can  drastically  change  the  thermodynamic  structure  of  the 
STBL  and  lead  to  evaporation  of  the  cloud  (Moeng  2000).  It  is  defined  as  the  process  that 
fluid  is  exchanged  across  a  density  interface  bounding  a  region  of  turbulent  flow  (Nichols 
and  Turton,  1986).  The  warm  and  dry  air  is  entrained  into  the  cloud  layer,  forming  an 
entrainment  parcel.  The  thermodynamics  of  the  entrainment  parcel  is  quite  complicated 
and  plays  a  key  role  in  STBL  evolution.  Mixing  with  the  dry  inversion  air  results  in 
sub  saturation  of  the  entrainment  parcel  and  thus  the  evaporation  of  the  cloud  droplets. 
As  a  result,  depletion  of  the  cloud  water  was  observed  in  the  entrainment  parcel  (Wang 
and  Albrecht  1994).  When  conditions  are  right,  the  entrainment  parcels  become 
negatively  buoyant  due  to  evaporative  cooling,  contributing  to  positive  buoyancy  flux  in 
the  entrainment  parcel.  This  process  is  referred  to  as  buoyancy  reversal.  Buoyancy 
reversal  strengthens  the  entrainment  by  driving  small  scale  turbulence  near  the  cloud  top, 
leading  to  positive  feedback  of  entrainment,  called  cloud-top  entrainment  instability 
(CTEI).  It  was  hypothesized  that  CTEI  may  cause  rapid  breakup  of  the  cloud  deck,  but 
such  rapid  breakup  has  not  been  verified  from  observations.  However,  Moeng  (2000) 
suggests  that  buoyancy  reversal  dominantly  controls  the  cloud  amount. 

The  entrainment  rate  is  considered  the  most  important  parameter  to 
determine  the  STBL’s  thermodynamic  and  cloud  structure  in  climate  models. 
Unfortunately,  this  variable  is  extremely  difficult  to  quantify.  The  Dynamics  and 
Chemistry  of  Marine  Stratocumulus  (DYCOMS-II)  represents  the  most  recent  organized 
effort  to  quantify  the  entrainment  rate.  The  reader  is  referred  to  Stevens  et  al.  (2003)  and 
the  many  DYCOMS-II  related  articles  on  the  latest  development  regarding  entrainment 
rate. 


e.  Wind  and  Surface  Fluxes 

Although  stratocumulus  clouds  are  often  found  over  the  cool  upwelling 


coastal  water,  surface  buoyancy  flux  maybe  positive  depending  on  the  air-sea 
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temperature  difference  and  thus  contribute  to  turbulent  mixing  in  the  STBL  (Wang  1993, 
de  Roode  and  Duynkerke  (1996).  In  general,  the  surface  buoyancy  flux  is  less  than 
10  Wm'2,  if  not  slightly  negative.  Over  relatively  warm  ocean,  the  surface  buoyancy  flux 
is  the  source  of  turbulence  in  the  decoupled  STBL  (de  Roode  and  Duynkerke,  1996) . 

The  surface  moisture  flux  is  the  main  source  of  moisture  into  the  BL.  However, 
compared  to  the  buoyancy  reversal  process,  the  net  impact  on  cloud  fraction  and  liquid 
water  path  from  the  surface  is  insignificant  (Moeng  2000).  The  surface  latent  heat  flux 
(LHF)  is  defined  as  the  flux  of  heat  from  water  evaporating  at  the  Earth’s  surface  and 
condensing  as  water  vapor  in  the  BL,  with  constant  temperature.  LHF  is  determined  by 
the  surface  relative  humidity,  temperature,  and  wind  speed  (Wood  2012).  In  contrast  to 
LHF,  the  surface  sensible  heat  flux  (SHF)  is  defined  as  the  conductive  heat  flux  from  the 
Earth’s  surface  to  the  atmosphere  and  is  dependent  on  a  change  of  temperature.  Both 
LHF  and  SHF  are  a  part  of  many  key  processes  transporting  energy  from  the  Earth’s 
surface  to  the  atmosphere.  All  physical  process  in  the  STBL  may  affect  the  various 
surface  fluxes  and  in  turn  affect  the  source  of  moisture  in  the  BL. 

Vertical  shear  of  the  mean  wind  in  STBL  are  observed  near  the  surface 
and  often  in  the  inversion  atop  the  cloud  layer.  Brost  et  al.  (1982)  presented  cases 
measured  from  the  California  coast  where  the  driving  force  of  turbulence  are  in  boundary 
layer  wind  shear,  not  radiative  cooling  in  the  STBL.  The  shear  in  the  inversion  may 
modify  entrainment  mixing  at  the  cloud  top.  Wang  et  al.  (2008)  and  Moeng  et  al.  (2005) 
both  examined  the  shear  effects  across  the  STBL  top  using  large  eddy  simulations  (LES). 
They  found  that  the  local  wind  shear  from  individual  eddies  can  also  modify  the 
entrainment  process  even  in  the  case  when  the  mean  wind  shear  is  negligibly  small. 

f  Precipitation 

Frequently,  stratocumulus  clouds  produce  light  precipitation  in  the  form  of 
drizzle  (Nicholls  and  Leighton  1986;  Wood  2012).  Precipitation  affects  cloud  fraction  by 
removing  water  from  the  BL  as  drizzle  at  the  surface  and  by  moistening  and  cooling 
below  the  cloud  as  the  falling  drops  evaporate  (Pincus  et  al.  1997).  Drizzle  is  found  20- 
40%  of  the  time  in  persistent  marine  stratocumulus.  VanZanten  et  al.  (2004)  suggests  that 
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drizzle  is  more  prevalent  than  previously  thought  and  can  even  induce  a  transition  in 
cloud  structure.  Higher  precipitation  rates  are  not  due  to  a  change  in  shape  of  the  drizzle 
droplets,  but  are  mainly  caused  by  more  frequent  occurrence  of  larger  drizzle  drops.  The 
rate  of  precipitation  is  higher  at  the  cloud  base  than  at  the  top,  so  it  is  often  characterized 
by  the  cloud  base  drizzle  rate.  Drizzle  stabilizes  the  STBL  by  warming  the  cloud  layer. 
This  also  reduces  turbulent  mixing  and  increases  stratification.  Below  the  cloud  base, 
drizzle  quickly  evaporates  due  to  the  small  droplet  size.  The  evaporation  of  drizzle  in  the 
sub-cloud  layer  is  rather  significant,  even  for  shallow  BL.  Over  both  the  land  and  the 
ocean,  the  maximum  amount  of  drizzle  occurs  in  the  early  morning  when  the  cloud 
thickness  is  greatest  (Wood  2012).  Sharon  et  al.  (2006)  observed  substantial  differences 
in  drizzle  characteristics  between  the  solid  stratocumulus  deck  and  rift  areas.  Drizzle  is 
clearly  more  evident  from  observations  in  rifts  than  the  surrounding  solid  cloud. 

3.  Diurnal  Variability 

Stratocumulus  clouds  exhibit  a  diurnal  pattern,  controlled  by  variation  in  solar 
insolation  throughout  the  day.  Maximum  cloud  coverage  normally  occurs  at  night 
(Skupniewicz  et  al.  1991).  Betts  (1990)  examined  the  diurnal  variation  in  the  BL 
thermodynamic  structure  off  the  CA  coast.  In  the  morning  after  sunrise,  the  solar  heating 
of  the  cloud  layer  results  in  decoupling  and  raises  the  cloud  base.  The  cloud  top  does  not 
rise  significantly,  the  cloud  is  thinned  and  may  result  in  fractional  cloudiness  (Pincus  et 
al.  1997).  Thus,  the  cloud  is  thinnest  when  the  maximum  amount  of  solar  energy  is 
absorbed.  When  the  cloud  is  thinnest,  it  transmits  the  most  solar  radiation  to  the  surface 
(Hanson  and  Gruber  1982). 
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III.  UPPEF  2012  MEASUREMENTS 


A.  UPPEF  2012  OVERVIEW 

UPPEF  2012  deployed  the  CIRPAS  Twin  Otter  in  stratocumulus  regimes  off  the 
CA  coast  in  order  to  examine  various  aspects  of  BL  cloud  and  turbulence 
parameterizations  for  applications  in  extended  range  forecast  models.  The  statistical 
parameterization  method  using  probability  distribution  functions  (PDF)  of  dynamic  and 
thermodynamic  variables  was  targeted  for  cloudy  boundary  layers  with  various  cloud 
fractions,  focusing  on  broken  stratocumulus  cloud  regimes.  Twelve  research  flights  by 
the  CIRPAS  Twin  Otter,  a  total  of  fifty  flight  hours,  were  dedicated  to  make  daytime 
comprehensive  aircraft  measurements  of  cloud  and  turbulence  in,  below,  and  above  the 
stratocumulus  cloud  off  the  CA  coast.  The  first  ten  flights  were  flown  for  the  cloud 
mission  and  the  last  two  flights  for  the  LIDAR  mission  west  of  Monterey  Bay  in  northern 
CA.  The  Controlled  Towed  Vehicle  (CTV)  was  utilized  for  five  of  the  twelve  flights 
(Table  1),  which  provided  a  new  measurement  capability  from  past  experiments.  The 
CTV  is  designed  to  characterize  surface  fluxes  as  low  as  10m  above  the  surface  while  the 
Twin  Otter  made  similar  measurements  at  its  flight  level.  However,  the  CTV  data 
collected  was  not  available  for  analysis  for  this  thesis  research. 

The  offshore  region  of  northern  and  central  CA  coast  is  one  of  the  most  extensive 
STBL  in  the  world,  making  it  an  excellent  location  to  study  stratocumulus  clouds. 
Focusing  on  the  multi-scale  variability  in  the  daytime  stratocumulus  cloud  provided  an 
opportunity  to  examine  the  relationship  between  multiple  processes  that  need  to  be 
incorporated  into  the  forecast  models.  Also,  since  the  clouds  tend  to  decouple  from  the 
surface  layer  during  the  afternoon,  based  on  their  diurnal  variability,  the  optimal 
measurement  time  for  studying  cloud  evolution  is  in  the  early  afternoon  along  the  CA 
coast.  Also,  the  CIRPAS  Twin  Otter  research  aircraft  is  based  at  Marina,  CA  (just  north 
of  Monterey,  CA)  and  has  easy  flight  range  of  the  stratocumulus  clouds.  A  graphical 
representation  of  the  flight  tracks  and  location  of  the  first  ten  flights  for  the  cloud  mission 
is  shown  in  Figure  9.  All  twelve  flights  are  summarized  in  Table  1,  to  include  the  flight 
number,  date,  time,  duration,  flight  pattern,  and  any  additional  comments.  A  variety  of 
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flight  patterns,  such  as  the  L-pattem,  U-pattem,  S-pattem,  box,  and  long  linear  transects, 
were  used  during  UPPEF  2012  to  sample  the  different  boundary  layer  and  cloud 
characteristics.  Vertical  stacks,  spiral  and  slant  path  soundings,  and  zig  zag 
horizontal/vertical  transects  were  also  utilized  in  flight  to  sample  mesoscale  and  small 
scale  variability  within  the  boundary  layer,  mean  radiation  measurements  above  and  in 
cloud,  low  level  turbulence,  and  microphysics  within  the  cloud  layer. 


UPPEF2012  Twin  Otter,  All  Cloud  Study  Flights 


Figure  9.  Flight  tracks  of  the  first  ten  UPPEF  flights  (RF01-RF10). 


Each  UPPEF  flight  was  airborne  for  approximately  3.5  to  6  hours  (see  Table  1  for 
exact  duration).  Using  the  flight  GPS  measurements,  individual  legs  are  defined  for  full 
soundings,  short  soundings,  and  level  legs.  Full  soundings  are  defined  as  legs  from  the 
surface  or  near  surface  to  above  the  cloud  top  (and  vice  versa).  Short  soundings  are 
defined  as  legs  originating  above  the  boundary  layer  top  going  into  the  boundary  layer 
(vice  versa)  and  level  legs  were  the  transits  with  the  same  altitude  at  constant  heading 
(normally  for  five  minutes  or  longer).  Figure  11a  shows  the  time  series  of  aircraft  altitude 
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(m)  and  heading  (°)  during  RF01  overlaid  each  full  sounding  during  flight.  The  full 
soundings  are  numbered  consecutively,  by  order  of  time  sequence.  Their  surface 
projections  are  shown  in  color  lines  overlaid  on  the  entire  flight  track  (Figure  1  lb).  In  the 
example  shown  in  Figure  1 1  (RF01),  five  full  soundings  were  completed.  Figure  12  is  the 
same  as  Figure  11,  except  it  is  for  the  level  legs  in  RF01.  Similar  flight  information  was 
created  and  different  types  of  legs  were  defined  for  all  of  the  other  eleven  UPPEF  2012 
flights. 


UPPEF12  Full  Soundings  RFQ1  20120831 


UPPEF12  Full  Soundings  RF01  20120831 


Figure  10.  (a)  Time  series  of  height  (m)  and  heading  (°)  showing  full  soundings  overlaid 

on  flight  track;  (b)  Latitude/longitude  plot  showing  full  soundings  overlaid  on 

flight  track  for  RF01 . 


UPPEF12  LL12  z  =  22  RF01  20120831 


UPPEF  12  Level  Legs  RF01  20120031 


Figure  11.  (a)  Time  series  of  height  (m)  and  heading  (°)  showing  level  legs  overlaid  on 

flight  track;  (b)  Latitude/longitude  plot  showing  level  legs  overlaid  on  flight 

track  for  RF01. 
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Table  1.  UPPEF  2012  Flight  Summary 


Flight  # 

Date 

Takeoff- 
Landing  Time 
(PDT) 

Flight 

Duration 

Major  Flight 
Pattern 

Comments 

RF01 

8/31/2012 

1207-1705 

5 

L-pattem 

Sampled  at  the  clear/c  loud  boundary,  ~  20-30  kts 
wind.  No  CTV. 

RF02 

9/4/2012 

1210-1720 

5.2 

S-pattem 

Mesoscale  variability  in  solid  cloud  cover  at  the 
beginning.  Broken  cloud  layer.  Multiple  sounding 
were  made.  ~10  kts  wind.  No  CTV. 

RF03 

9/6/2012 

1228-1717 

4.8 

S-pattem 

Thin  solid/broken  cloud  cover.  Weak  and 
sometimes  no  wind,  high  SST.  Twin  Otter  had  a 
hard  time  staying  in  cloud  because  the  cloud  layer 
was  very  thin.  No  CTV. 

RF04 

9/7/2012 

0941-1500 

5.3 

U-pattern 

Gradual  increase  of  cloud  fraction  from  the  coast 
to  the  west  side  of  the  U-pattern.  U-pattem  at 
100  ft  above  surface  for  turbulence  and  cloud 
survey.  Vertical  stack  in  the  cross-win  direction 
at  four  levels  below  the  cloud  top.  More  clearing 
later  in  the  flight.  No  CTV. 

RF05 

9/11/2012 

0953-1310 

5.3 

U/box-pattem 

Solid  cloud  with  visible  roll  structure.  Large 
clearing  to  the  north.  Some  strong  wind  to  the 
west  in  the  alongwind  leg.  Variable  surface 
waves  and  swell.  No  CTV. 

RF06 

9/13/2012 

1137-1507 

3.5 

Long  linear 

transect 

1-hour  long  straight  leg  to  the  west  to  capture 
mesoscale  features  and  broken  clouds  further 

offshore.  First  mission  of  UPPEF  with  CTV. 

RF07 

9/17/2012 

1022-1537 

5.3 

U-pattern 

Persistent  solid  cloud  with  clear  mesoscale 
features.  Cross-wind  leg  of  U-pattern  is  90k  to 
capture  mesoscale  features.  No  CTV. 

RF08 

9/19/2012 

1102-1617 

5.3 

Long  linear 

transect 

Long  linear  transect  going  through  cloud  (near 
coast)  to  clear  transition.  Roll  structure.  CTV 
used. 

RF09 

9/20/2012 

0910-1430 

5.3 

U-pattern 

Thin  cloud,  dissipated  later  on.  Most  of  the 
measurements  were  made  in  clear  with  some  thin 
puffy  clouds.  CTV  used. 

RF10 

9/26/2012 

1006-1342 

3.6 

Linear 

UPPEF  cloud  mission  with  CTV  and  with 
TODWL  installed.  Maximum  flight  duration  is 
3.6  hr  in  this  configuration.  Solid  to  broken  cloud 
transition  from  coast  to  further  west. 

RF11/RF12 

9/30/2012 

0920-1636 

5.8 

Varied 

Two  flights  on  this  day  after  several  days  of 
waiting  for  clear  condition.  These  flights  are  not 
part  of  UPPEF  2012  cloud  missions,  but  are  part 
of  EDMF  study.  Both  with  CTV.  Wind  lidar 
data  not  available  for  most  of  second  flight  due  to 
high  cabin  temperature. 

Upper  air  soundings  are  important  for  future  modeling  efforts  to  characterize  the 
temperature,  humidity,  pressure,  and  wind  at  above  boundary  layer  levels.  For  this 
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purpose,  two  VIASALA  rawinsondes  were  launched  for  all  flights  from  CIRPAS  at  the 
Marina  municipal  airport,  except  for  the  first  two  flights  where  only  one  rawinsonde  was 
available.  Examples  of  temperature  and  dew  point  vertical  profiles  from  a  sounding  are 
shown  in  Figure  13  for  RF01,  RF02,  and  RF03.  Figure  14  displays  the  corresponding 
wind  speed  and  wind  direction  profiles.  The  inversion  base  and  any  significant  features 
are  annotated  in  each  plot.  As  seen  in  Figure  14a-d,  the  air  is  much  drier  above  the  BL 
than  below  as  indicated  by  the  increased  separation  between  temperature  and  dew  point. 
This  is  the  effect  of  the  ocean  surface  as  a  general  source  of  water  vapor  and  the  capping 
inversion.  The  rawinsonde  from  RF03  shows  the  air  column  from  4000-5000m  to  be 
nearly  saturated  (Figure  13c,  13d).  For  RF01,  the  temperature  inversion  was  at  768m  over 
Marina  (at  the  ascent)  and  808m  over  Salinas  Valley  (at  the  descent)  (Figure  13a).  The 
wind  was  northwesterly  at  the  surface  and  a  westerly  wind  maximum  (7  ms"1)  was  at 
200m  ASL.  The  wind  speed  varied  from  2-7  ms"1  for  the  entire  rawinsonde  launch  and 
the  wind  direction  shifted  from  westerly  in  the  BL  to  southerly  above  (Figure  14a).  The 
temperature  inversion  for  RF02  was  at  360m,  half  the  altitude  of  the  first  flight  (Figure 
13b).  The  launches  from  RF03,  suggest  high  temporal  variability  of  the  BLH.  The  height 
of  the  inversion  for  the  first  launch  was  at  429m  (Figure  13c)  while  the  height  4.5  hours 
later  was  suppressed  to  250-284m  (Figure  13d). 

There  were  two  types  of  launches  for  UPPEF  2012.  The  first  was  designed  to  go 
‘up  and  down’  in  order  to  get  a  good  depiction  of  the  boundary  layer  and  the  downward 
soundings  provide  the  low-level  measurements  away  from  any  ship  effects.  Examples  of 
this  type  are  the  balloon  launch  for  RF01  and  the  second  launch  during  RF03,  as 
indicated  in  Figure  14.  The  second  type  is  a  synoptic  sounding,  to  go  “up  only,”  in  order 
to  get  data  from  the  entire  troposphere.  Examples  are  the  balloon  launches  for  RF02  and 
the  first  for  RF03.  In  the  wind  profiles,  the  upper  level  jet  is  seen  in  the  synoptic 
soundings  (Figure  14b,  14c),  which  the  “up  and  down”  profiles  for  low  levels  did  not 
capture  (Figure  14a,  14d).  The  mid-level  jet  on  the  two  days  has  similar  magnitude, 
although  the  vertical  variation  leading  to  the  jet  core  is  rather  different  on  these  two  days. 
Even  though,  the  type  of  sounding  is  labeled  “up  only”  data  was  still  recorded  as  the 
rawinsonde  descended  in  the  atmosphere  once  the  balloon  was  popped  (no  parachutes 
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were  used  for  the  descent).  Similar  plots  were  created  for  every  rawinsonde  launch  during 
each  flight. 


RF01  Rawinsonde  Launch  2103Z  RF02  Rawinsonde  Launch  1700Z 


RF03  Rawinsonde  Launch  1500Z 


RF03  Rawinsonde  Launch  1941Z 


Figure  12.  Annotated  T  and  Td  profiles  for  (a)  RF01,  (b)  RF02,  (c)  RF03,  and 
(d)  the  second  sounding  of  RF03  from  rawinsonde  launches. 
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RF01  Rawinsonde  Launch  2103Z  RF01  Rawinsonde  Launch  2103Z  RF02  Rawinsonde  Launch  1700Z  RF02  Rawinsonde  Launch  1700Z 


RF03  Rawinsonde  Launch  1500Z  RF03  Rawinsonde  Launch  1500Z  Rp°3  Rawinsonde  Launch  1941Z  RF03  Rawinsonde  Launch  1941Z 


Figure  13.  Annotated  wind  speed  and  direction  vertical  profdes  for  (a)  RF01, 

(b)  RF02,  (c)  RF03,  and  (d)  the  sounding  of  RF03  from  rawinsonde  launches. 


B.  CIRPAS  TWIN  OTTER  INSTRUMENTATION  AND  MEASUREMENTS 

The  CIRPAS  Twin  Otter  was  outfitted  with  improved  instrumentation  for  UPPEF 

2012.  There  were  four  different  GPS  receivers,  owned  by  CIRPAS  and  UCI,  that 

provided  flight  path  and  navigation  data.  An  improved  5-hole  Radome  gust  probe 

measured  wind  and  turbulence  was  installed.  There  were  two  fast-response  Rosemont 

total  temperature  sensors  from  UCI  and  CIRPAS  for  redundancy  and  accuracy.  Other 

sensors  measured  static  and  dynamic  pressure,  dew  point,  water  vapor,  absolute 

humidity,  and  carbon  dioxide  concentration.  There  were  various  upward  and  downward 

looking  pyranometers  measuring  SST,  up-  and  down-welling  solar  irradiance,  up-  and 

down-welling  IR  irradiance,  and  IR  temperature.  There  were  also  sensors  collecting 

aerosol  and  cloud  properties,  such  as  the  Passive  Cavity  Aerosol  Spectrometer  Probe 

(PCASP)  and  the  Forward  Scatter  Spectrometer  Probe  (FSSP).  The  measurement  was 
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monitored  in  real-time  at  CIRPAS  through  satellite  links.  The  meteorological  data  was 
processed  and  made  available  for  the  user  from  CIRPAS  at  1  Hz  and  10  Hz.  Dr.  Djamal 
Khelif,  UC  Irvine,  also  logged  some  of  the  fast  sampling  data  separately  and  made  them 
available  for  this  research  at  40  Hz.  A  photograph  of  the  front  of  the  Twin  Otter  is  shown 
in  Figure  15  annotated  with  the  locations  of  the  major  instruments  used  during  UPPEF 
2012.  The  CIRPAS  and  UCI  instruments  used  for  UPPEF  2012,  along  with  their 
description,  measured  variable(s),  and  data  accuracy,  are  listed  in  Table  2.  The  onboard 
instrumentation  is  similar  to  previous  Twin  Otter  research  flights,  as  described  by  Khelif 
et  al.  (1999)  and  Khelif  et  al.  (2005). 


Figure  14.  A  photograph  of  CIRPAS  Twin  Otter  showing  placement  of  major 
instruments  used  during  UPPEF  2012  (Courtesy  of  Djamal  Khelif, 
University  of  California,  Irvine) 
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Table  2.  CIRPAS  Twin  Otter  Instrument  List  used  for  UPPEF  2012  (input  provided 
by  Haflidi  Jonsson,  Djamal  Khelif,  and  Anthony  Bucholtz). 


Instrument 

Purpose 

Description 

Range 

TANS  Vector  platform  attitude 
(CIRPAS) 

Flight  Path/  Navigation 

GPS,  Pitch,  Roll,  Heading 

NovAtel  GPS  (CIRPAS) 

Fhght  Path/  Navigation 

GPS,  Lat,  Lon,  Alt,  ground  speed 
and  track 

C-MIGITS-III,  GPS/INS  (2  units, 
CIRPAS/UCI) 

Flight  Path/  Navigation 

Lat,  Lon,  Alt,  ground  speed  and 
track,  pitch,  roll,  heading 

OXTS-RT3003  GPS/INS  (UCI) 

Fhght  Path/  Navigation 

Lat,  Lon,  Alt,  ground  speed  and 
track,  pitch,  roll,  heading 

had  a  OXTS  base  station  to  correct 
for  post-flight  differential  GPS 

Rosemount  Total  Temperature 
sensor  (CIRPAS) 

Meteorological  parameters 

Total  Temperature 

-50  to  +50  °C 

Rosemount  Total  Temperature 
sensor  (UCI) 

Meteorological  parameters 

Total  Temperature 

-60  to +40  °C 

Edgetech  Vigilant  Dew  Point 
Temperature  sensor  (CIRPAS) 

Meteorological  parameters 

Chilled  Mirror  device 

-50  to  +50  °C 

Setra  270  transducer-  Static 
pressure  sensor  (CIRPAS) 

Meteorological  parameters 

Static  (barometric)  pressure 

600-1 100  mb +/-  .25  mb 

Setra  270  transducer  -  total  pressure 
sensor  (CIRPAS) 

Meteorological  parameters 

total  (static  +  dynamic)  pressure 

600-1 100  mb  +/-  .25  mb 

Wind  -  Radome,  flow  angle  probe 
(CIRPAS) 

Meteorological  parameters 

Angle  of  attack  and  angle  of 
sideshp  differential  pressures 

-37.3-37.3  mb  +/-  .1  mb 

Heiman,  Heitronics  KT  19.85 
Pyrometer  (CIRPAS) 

Meteorological  parameters 

Sea  surface  temperature 

-5  to  45  °C  adjustable 

Heitronics  KT19.85  Upward-looking 
IRTemp  (UCI) 

Meteorological  parameters 

IR  Temperature 

-50  to +20  °C 

Up-looking  modified  Kipp  &  Zonen 
CM-22  pyranometer  (CIRPAS/NRL) 

Meteorological  parameters 

Down-welling  Solar  Irradiance 

0-1400  Wm2 

Down-looking  modified  Kipp  & 
Zonen  CM-22  pyranometer 
(CIRPAS/NRL) 

Meteorological  parameters 

Up-welling  Solar  Irradiance 

0-1400  Wm'2 

Up-looking  modified  Kipp  &  Zonen 
CG-4  pyranometer  (CIRPAS/NRL) 

Meteorological  parameters 

Down-welling  IR  Irradiance 

0-600  Wm'2 

Down-looking  modified  Kipp  & 
Zonen  CG-4  pyranometer 
(CIRPAS/NRL) 

Meteorological  parameters 

Up-welling  IR  Irradiance 

0-600  Wm2 

Fast  LI-COR  7500  H20/C02 
Analyzer  (UCI) 

Meteorological  parameters 

Water  vapor  and  C02  densities 

adjustable 

Fast  UCI-  Modified  KH20  (UCI) 

Meteorological  parameters 

Absolute  Humidity 

adjustable 

Passive  Cavity  Aerosol 
Spectrometer  Probe  (PCASP)  - 
SPP200  (CIRPAS) 

Aerosol  properties 

Particle  Size  Distribution 

0.1  -3.0  pm  (100-3000  nm) 

2  Condensation  Particle  Counters 
(TSI 3010)  (CIRPAS) 

Aerosol  properties 

Particle  concentration 

Dp  >  0.01  pm  (10  nm);  Dp  >  0.015  pm 
(15  nm) 

Ultrafine  Condensation  Particle 
Counter  (TSI  3025)  (CIRPAS) 

Aerosol  properties 

Particle  concentration 

Dp  >0.003  pm  (3  nm) 

Forward  Scatter  Spectrometer  Probe 
(FSSP)  -  SPP100  (CIRPAS) 

Cloud  properties 

Particle  Size  Distribution 

2.0  -  30.0  pm 

Cloud,  Aerosol,  and  Precipitation 
Spectrometer  (CAPS)  (CIRPAS) 

Cloud  properties 

Particle  and  Drop  Size 
Distribution 

0.5-50  pm  (scatter) 

25-1600  pm 

PVM-100A  probe  (CIRPAS) 

Cloud  properties 

Liquid  Water  Content  of  clouds 
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C.  SYNOPTIC  CONDITIONS  DURING  UPPEF  2012 

The  synoptic  conditions  during  the  UPPEF  2012  measurement  period  are 
summarized  in  this  section.  On  August  31,  2012,  the  synoptic  pattern  along  the  West 
Coast  featured  a  quasi-stationary  500-mb  trough  axis  situated  along  the  coast  in  a  north- 
south  orientation,  causing  strong  onshore  flow.  The  EPAC  High  was  located  to  the  west 
of  California  and  Oregon  border,  and  moving  slowly  ENE.  During  the  night  before,  as 
the  atmosphere  stabilized,  extensive  solid  marine  stratocumulus  with  small  horizontal  cell 
structures  developed  all  along  the  coast  from  Cape  Mendocino,  south  and  coastward  of 
125W.  West  of  125W,  there  were  very  deep  stratocumulus  clouds  with  large  cell 
structures  with  a  narrow  clear  gap  separating  the  two  marine  stratocumulus  regimes. 
Along  the  coast,  cell  sizes  were  on  the  order  of  a  few  kilometers,  while  west  of  125W  the 
cell  sizes  were  around  50km  in  diameter.  During  RF01,  the  wind  was  primarily 
northwesterly  in  the  BL.  The  Ft.  Ord  profiler  showed  a  deep  marine  layer,  with  the 
inversion  near  762m  (2500ft). 

Over  the  next  few  days,  the  upper  level  trough  shifted  its  orientation  to  northeast- 
southwest.  The  EPAC  High  built  to  1032mb,  moved  north  to  the  west  of  the  Oregon  and 
Washington  border,  and  tightened  the  surface  pressure  gradient  over  northern  CA  coast. 
A  thermal  trough  set  up  over  Nevada  and  northern  CA  extending  from  a  thermal  low  over 
New  Mexico.  Tropical  Storm  (TS)  John  was  a  weak  system  that  formed  off  the 
southwestern  coast  of  Mexico  on  September  2,  2012.  TS  John  weakened  over  the  cool 
waters  west  of  the  Baja  Peninsula  on  4  September  to  a  Tropical  Depression  (TD)  and 
advected  midlevel  moisture  into  the  UPPEF  region.  For  the  second  flight,  on  September 
4,  2012,  a  weak  mid-  to  upper  level  low  remained  west  of  the  northern  CA  coast,  with  an 
upper  level  ridge  over  the  Rockies.  There  was  a  stable  marine  layer  with  heights  over 
most  of  the  coast  near  365m  (1200ft)  and  cloud  cover  extended  well  west  of  Monterey 
Bay  past  125W.  The  cloud  bases  were  quite  low  at  sunrise  and  rapidly  rose  during  the 
morning  hours.  By  noon,  clearing  was  evident  everywhere  inland  of  the  beach  except  for 
the  first  couple  of  miles  inland  around  Monterey  and  Marina.  The  surface  winds  were 
northwesterly  at  5-7ms_1  offshore  of  Monterey  Bay,  while  above  the  BL  the  winds  shifted 
to  southwesterly. 


32 


On  5-6  September  2012,  a  cut  off  upper  level  low  formed  west  of  Monterey  Bay 
and  moved  slowly  northeast.  This  movement  pushed  the  weakened  EPAC  high  north  into 
the  Gulf  of  Alaska.  For  the  third  flight,  on  September  6,  2012,  the  stable  marine  layer 
along  the  central  CA  coast  deepened  with  an  inversion  height  approximately  430m.  By 
midday,  the  cloud  layer  eroded  out  from  the  coastline  to  nearly  50nm  offshore.  The  cloud 
bases  were  initially  near  60m  (200ft)  and  lifted  as  the  clouds  in  the  nearshore  region 
burned  off.  The  western  edge  of  a  midlevel  tropical  moisture  plume  was  located  just  west 
of  Monterey  Bay  (also  seen  in  the  rawinsonde  plots  in  Figure  13).  Between  the  coastline 
and  130W,  a  weak  surface  eddy  or  gyre  developed,  with  northeasterly  surface  wind  flow 
from  the  Washington  and  Oregon  border  down  to  Cape  Mendocino.  South  of  Cape 
Mendocino,  off  Monterey  Bay,  the  wind  was  west-northwesterly  in  the  BL  and  southerly 
above.  There  was  a  weak  cyclonic  structure  in  the  stratocumulus  field  from  35N  and 
north  to  the  Oregon  border. 

On  7  September  2012,  for  the  fourth  flight,  the  surface  pressure  gradients  along 
the  coast  relaxed  to  a  weak  onshore  pressure  gradient  with  little  north-to-south  gradient. 
Winds  were  northwesterly  in  the  BL,  while  above  there  was  a  weak  southerly  flow 
lowering  the  inversion  height  to  around  250m.  In  the  morning,  there  was  strong  surface 
warming  due  to  clear,  sunny  skies  along  the  coast  extending  westward  to  the  ragged  edge 
of  the  solid  stratocumulus  field  roughly  along  125W.  The  cloud  boundary  remained 
quasi-stationary  throughout  the  night  and  appeared  to  be  retreating  slowly  in  the  early 
morning  hours.  By  8  September  2012,  the  EPAC  high  weakened  to  1020mb  and  then 
reformed  west  of  southern  CA.  The  following  day,  an  amplified  upper  level  ridge  set  up 
over  northern  Montana  with  a  deep  trough  west  of  California.  The  longwave  pattern  then 
moved  east  and  both  the  ridge  and  trough  decreased  in  amplitude.  By  the  afternoon  of 
1 1  September,  and  on  the  fifth  flight,  the  EPAC  high  built  to  1031  mb  and  moved  north  of 
CA  toward  the  Gulf  of  Alaska.  A  thermal  trough  set  up  over  northern  CA  which  extended 
from  the  thermal  low  over  Arizona  (desert  region).  The  large  surface  EPAC  high  pressure 
in  the  Gulf  of  Alaska  brought  strong  northeasterly  winds  along  the  coast  southward  from 
the  WA/OR  border  to  just  north  of  San  Francisco  Bay.  A  wind  reversal  event  with  weak 
cyclonic,  coastally  trapped  southerly  flow  occurred  right  along  the  immediate  coast  of 
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Monterey  County.  Further  offshore,  the  winds  were  northwesterly  in  the  BL  along 
RF05’s  flight  track.  The  inversion  rose  from  the  previous  flight  to  400-500m  along  the 
coast. 

The  upper  level  ridge  moved  into  the  Gulf  of  Alaska  by  12  September  0600  UTC 
and  continued  to  move  east.  A  thermal  trough  set  up  again  over  northern  CA  by 
1800UTC.  On  September  13,  2012,  the  sixth  flight,  a  strong  500-mb  ridge  has  amplified 
along  the  West  Coast.  The  wind  reversal  continued  with  the  southerly  surge  of  cool, 
dense  marine  air  that  propagated  northward  out  of  the  Monterey  area.  The  Monterey  area 
was  under  an  area  of  strong,  warm  subsidence  which  has  kept  the  marine  layer  along  the 
coast  compressed  to  approximately  300-350m.  The  winds  were  southeasterly  above  the 
inversion  and  light  west-northwesterly  in  the  MBL.  There  was  an  extensive 
stratocumulus  deck,  evident  in  satellite  imagery,  continued  to  hug  the  coastline.  The  early 
morning  cloud  bases  were  very  low  with  dense  fog  inland.  By  early  afternoon,  the  fog 
lifted  and  burned  back  toward  the  beaches. 

During  the  next  few  days,  the  elongated  EPAC  High  (1028mb)  moved  south  from 
the  Gulf  of  Alaska  toward  the  Washington  and  Oregon  border.  The  upper  level  ridge 
moved  eastward  away  from  the  CA  coast.  On  September  15,  a  weak  trough  over  the 
EPAC  moved  east  toward  the  west  coast.  By  September  16,  a  300mb  ridge  was  over 
San  Francisco  and  Monterey  Bay  area.  During  the  seventh,  eighth,  and  ninth  flights, 
17-20  September  2012,  a  Rex  block  formed  with  an  upper  level  low  west  of  the  OR/C  A 
border  and  an  upper  level  ridge  to  the  north.  The  boundary  layer  initially  deepened  to 
about  600m  and  widespread  solid  and  broken  stratocumulus  coverage.  There  was  little 
change  in  cloud  cover  with  a  very  strong  inversion  that  kept  the  stratocumulus  layer 
firmly  capped  for  the  seventh  flight.  There  was  strong  southerly  flow  aloft  and  light 
surface  westerly  winds  along  the  coast.  During  the  eighth  flight,  the  solid  stratocumulus 
cloud  became  broken  and  burned  off  inland  during  the  afternoon,  with  the  temperature 
inversion  diffusing  and  falling  to  400m  along  the  coast.  The  marine  layer  was  not  well 
defined  for  this  flight.  The  upper  level  southerly  winds  weakened  substantially  from  the 
past  flight  and  surface  winds  shifted  northwesterly.  During  the  ninth  flight,  there  were 
mostly  clear  skies  across  the  region,  with  small  areas  of  stratus  over  coastal  waters,  with 
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similar  inversion  heights  of  350-400m.  Winds  were  northwesterly  in  the  BL  and  westerly 
above. 

By  September  23  at  0600  UTC,  the  upper  level  low  in  the  EPAC  pushed  through 
the  Rex  block  and  became  cut  off,  as  it  was  under  the  amplified  upper  level  ridge  and 
then  moved  slowly  southeast  over  the  next  few  days  until  it  merged  with  the  longwave 
pattern.  On  the  tenth  flight,  on  26  September  2012,  a  highly  amplified  ridge  covered  the 
eastern  Pacific  Ocean  with  a  deep  trough  over  eastern  North  America.  The  deepened 
marine  layer  (400-5 00m)  and  onshore  flow  allowed  stratus  to  penetrate  inland  and  along 
the  coast  for  the  duration  of  the  flight.  At  500-mb,  the  decaying  Rex-block  with  a  high 
centered  west  of  Baja  California  coupled  with  a  weak  low  centered  west  of  WA/OR 
border  displaced  the  EPAC  high  south  of  its  normal  position  and  induced  a  cyclonic 
circulation  in  the  stratocumulus  cloud  field.  For  the  last  two  flights,  on  30  September 
2012,  another  Rex  blocking  pattern  affected  the  Monterey  Bay  region,  which  was  under  a 
strong  500mb  ridge.  This  strong  ridge  compressed  the  inversion  along  with  warmer  air 
aloft  into  a  surface  based  inversion,  less  than  a  hundred  meters  deep  along  the  coast;  and 
even  washed  out  the  marine  layer  in  many  places,  allowing  for  clear  skies  in  the  area. 
Winds  were  westerly  at  the  surface  veering  to  easterly  aloft  (up  to  12000m). 

D.  CLOUD  STRUCTURE  VARIABILITY 

Throughout  September  2012,  the  cloud  pattern  varied  from  clear,  broken,  and 
solid  stratocumulus  cloud  depending  on  the  boundary  layer  and  synoptic  patterns  that 
occurred.  The  MODIS  images  are  shown  in  Figure  15  and  represent  the  cloud  conditions 
at  approximately  2145UTC  for  each  flight,  except  for  RF11.  RFll’s  flight  occurred 
earlier  in  the  morning  and  the  MODIS  time  did  not  cover  that  time,  so  the  GOES- 15 
visible  image  was  used  instead.  During  RF01,  there  was  extensive  stratocumulus  cloud  in 
the  morning  in  Monterey  Bay  and  along  the  coast.  Throughout  the  day,  the  cloud 
dissipated  from  the  Bay  and  to  the  south.  The  second  flight  had  mainly  solid 
stratocumulus  cloud,  with  some  breaks  in  cloud  during  the  end  of  the  flight.  The  clouds 
were  broken  in  the  S-pattem  for  the  third  flight,  and  clear  near  the  coast.  The  cloud 
pattern  for  the  fourth  flight  had  larger  cellular  clouds  in  the  westernmost  legs  and 
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primarily  clear  skies  near  the  coast.  The  fifth  flight  appeared  to  have  rifts  in  the  western 
part  of  the  track  with  solid  cloud  towards  the  coast.  A  long  linear  transect  was  flown  for 
the  sixth  flight  capturing  the  solid  cloud  near  the  coast,  broken  cloud  in  the  middle,  with 
more  solid  cloud  to  the  west.  The  seventh  flight  had  primarily  large  cellular  structures  in 
the  U-pattem  track,  clearing  near  the  coast.  There  was  a  strip  of  stratocumulus  cloud  in 
the  eight  flights,  with  clear  conditions  over  land  and  west  of  the  track.  The  ninth  flight 
had  fairly  clear  conditions  and  solid  stratocumulus  during  the  entire  tenth  flight.  The  last 
two  flights  were  mostly  clear,  with  few  thin  scattered  clouds  in  the  region. 


Figure  15.  MODIS  images  capturing  the  cloud  structure  variability  for  the  twelve 
flights  (RF01-RF12)  during  UPPEF  2012.  Blue  lines  indicate  flight  tracks. 

E.  MARINE  BOUNDARY  LAYER  HEIGHT 

According  to  Seidel  et  al.  (2010),  there  are  seven  different  methods  to  estimate 
PBL  heights:  the  parcel  method,  level  of  potential  temperature’s  maximum  vertical 
gradient,  the  base  of  the  elevated  temperatures  inversion,  the  top  of  the  surface  based 
inversion  (SBI),  the  level  of  the  minimum  vertical  gradient  of  specific  humidity  (q)  or 
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relative  humidity  (RH)  or  refractivity  (N).  The  parcel  method  evaluates  the  BL  height  by 
identifying  the  height  at  which  0V  is  the  same  value  aloft  as  at  the  surface  (Seibert  et  al. 
2000).  It  essentially  determines  the  level  where  a  hypothetical  parcel  of  air  will  be  in 
equilibrium  with  its  surrounding  environment.  The  base  of  an  elevated  temperature 
inversion  serves  as  a  cap  to  mixing  below.  The  top  of  a  SBI  indicates  a  stable  boundary 
layer  and  defines  the  top  of  the  PBL.  The  other  methods  assume  an  unstable  or  neutral 
BL  structure  and  are  not  used  for  stable  conditions  (Seidel  et  al.  2010).  The  method  used 
while  analyzing  the  UPPEF  2012  data  research  was  a  combination  of  all  the  methods,  but 
primarily  using  the  base  of  the  temperature  inversion,  cloud  top,  minimum  gradient  in  q, 
and  turbulence  measurements.  The  top  of  the  SBI  was  used  to  determine  the  BL  height 
for  the  last  two  flights.  A  combined  BL  height  graphic  of  each  UPPEF  2012  flight  is 
displayed  in  Figure  16.  Each  flight  is  represented  by  a  different  symbol.  The  varying 
heights  are  represented  by  different  colors,  as  shown  in  the  graph’s  color  bar.  Each  flight 
day  exhibits  a  downward  sloping  of  the  BL  height  toward  to  the  coast  from  the  open 
ocean,  as  described  in  the  conceptual  model  in  Chapter  II. 

RF01  had  the  highest  BL  heights  between  800-950m.  RF02’s  BL  heights 
decreased  by  about  half  and  started  near  350m  at  the  coast  increasing  to  525m  offshore. 
The  BL  heights  in  RF03  were  300m  along  the  coast  increasing  to  almost  700m  offshore. 
The  BL  heights  for  RF04  and  RF05  were  around  700m  offshore  and  at  the  coast 
approximately  400m.  During  RF06,  the  coastal  BLH  was  at  300m,  while  at  the  furthest 
west  leg  was  at  a  height  of  720m.  The  BL  heights  during  RF07  were  considerably  higher 
along  the  coast,  at  600m,  increasing  to  800m  offshore.  There  was  about  a  250m  change 
during  RF08  from  680m  offshore  near  124. 5W  to  430m  along  the  coast.  The  BL  heights 
near  123W  rose  considerably  from  RF08  from  600m  to  800m.  Along  the  coast  BLH  were 
similar  near  400m.  The  tenth  flight  had  a  BL  with  gradual  downward  slope  from  800- 
750m  from  1 24-1 25 W,  then  a  much  larger  slope  starting  at  750m  to  400m  to  the  coast. 
The  last  two  flights  RF1 1  and  RF12,  had  boundary  layer  heights  considerably  lower  than 
all  other  flights,  as  there  was  a  surface  based  inversion  on  this  cloud  free  day,  reducing 
the  heights  to  below  200m. 
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Figure  16.  Combined  BLH  for  UPPEF  2012  flights 
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IV.  BOUNDARY  LAYER  AND  CLOUD  CHARACTERISTICS 


A.  BOUNDARY  LAYER  HEIGHT  VARIABILITY 

A  common  pattern  in  each  flight  was  the  downward  slope  of  the  boundary  layer 
top  towards  the  east  (shoreward)  as  presented  in  Chapter  II,  Section  B.  The  height  of  the 
boundary  layer  top  changed  as  a  result  of  synoptic  and  mesoscale  pattern  changes.  Two 
of  the  ten  flights,  RF02  and  RF10,  were  good  examples  of  boundary  layer  variability  in 
solid  to  broken  cloud.  This  section  will  show  satellite  images  from  the  Geostationary 
Operational  Environmental  Satellite  (GOES)  West,  also  called  GOES- 15,  visible  sensor 
and  the  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  sensor  on  the  Terra 
polar  orbiting,  sun  synchronous  satellite,  with  resolutions  of  1km  and  250m,  respectively. 
Boundary  layer  heights,  vertical  profiles,  and  cross  sections  of  cloud  liquid  water  ( qc ), 
potential  temperature  ( 9 ),  specific  humidity  (q),  wind  speed,  and  wind  direction  from  full 
soundings  are  also  shown  in  this  section. 

1.  RF02  (04  September  2012) 

There  was  solid  to  broken  stratocumulus  cloud  during  the  second  flight,  RF02,  on 
4  September  2012.  Three  consecutive  GOES- 15  visible  satellite  images  off  Monterey 
Bay  (Figure  17),  showing  the  spatial  and  temporal  cloud  variability  for  the  duration  of  the 
flight.  The  cloud  remained  solid  to  broken  with  clearly  identifiable  mesoscale  structure 
especially  earlier  in  the  flight.  Total  flight  duration  was  5.2  hours  using  mainly  the  S- 
pattem  at  two  levels  in  cloud  and  two  levels  below.  A  zoomed  in  satellite  image  from  the 
MODIS  overlaid  on  Google  Earth  with  the  flight  track  and  color  coded  full  sounding  legs 
are  shown  in  Figure  18.  All  of  the  full  sounding  legs  went  through  cloud  top  and 
provided  valuable  data  in  evaluating  the  boundary  layer  structure.  Cloud  liquid  water 
(< qc ),  potential  temperature  (9),  specific  humidity  (q),  wind  speed,  and  wind  direction,  are 
plotted  for  each  full  sounding  of  RF02  in  Figure  19.  The  profiles  show  the  horizontal 
variability  in  the  boundary  layer.  The  cloud  top  is  lowest  near  the  coast  and  rises  offshore 
and  correlates  well  with  the  level  of  maximum  9  increase,  q  decrease,  and  sometimes  a 
change  in  wind  speed  or  direction. 
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Figure  17.  Three  sequential  GOES  15  satellite  images  on  4  September  2012  showing 
mostly  solid  cloud  during  RF02.  Blue  line  is  flight  track,  with  plane  icon 
indicating  location  of  TO  at  time  of  satellite  pass. 


Figure  18.  MODIS  image  on  4  September  2012  at  2045  UTC  overlaid  on  Google  Earth 
with  RF02  flight  track  and  the  horizontal  projection  of  all  full  soundings  (FS1- 
FS9,  colored  via  legend  in  Figure  19). 
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Figure  19.  Vertical  profiles  for  all  full  soundings  in  RF02,  showing  cloud  and  boundary 
layer  variability  in  qc,  6,  q,  wind  speed,  and  wind  direction.  The  spikes  in  FS  1 
and  FS9  are  due  to  interference  with  high  frequency  communications  onboard 
the  TO  during  takeoff  and  landing. 


The  boundary  layer  height  was  estimated  from  the  vertical  profiles  as  in  Figure 
19,  and  overlaid  on  the  flight  track  in  Figure  20.  In  the  beginning  of  the  flight,  the 
boundary  layer  height  was  269m  from  the  take-off  sounding.  The  height  progressively 
increased  offshore  to  526m  at  the  westernmost  part  of  the  track.  As  the  TO  headed  back 
to  shore  near  the  end  of  the  flight,  the  boundary  layer  decreased  towards  the  East.  The 
last  full  sounding  from  landing  was  at  375m,  which  is  higher  than  the  take  off  initial 
sounding,  which  is  likely  related  to  daytime  heating  of  the  surface  over  land. 
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Figure  20.  Boundary  Layer  heights  for  RF02.  Dotted  line  represents  flight  track. 

Colored  dots  and  numbers  represent  boundary  layer  height. 


Liquid  water  content,  specific  humidity,  temperature,  and  wind  speed  for  FS1- 
FS5  and  SSI  for  RF02  are  displayed  using  the  contour  plot  function  in  MATLAB,  in 
Figure  21.  The  red  colors  are  the  highest  values  and  blue  colors  the  lowest.  The  black 
solid  line  indicates  the  approximate  slope  for  the  BLH,  estimated  to  be  at  2  m  per 
kilometer.  The  figure  shows  the  upward  tilt  in  the  cloud  top  and  base  with  cloud 
thickness  slightly  increasing  offshore  and  thickest  near  122.9W  (FS4).  The  qc  increases 
from  the  base  of  the  cloud  with  maximum  values  (>0.4  gkg'1)  near  the  cloud  top.  Specific 
humidity  (Figure  21b)  is  high  and  fairly  uniform  with  height  in  the  boundary  layer  and 
through  the  cloud  top,  decreasing  rapidly  above  the  cloud  top  height.  Higher  moisture 
(7.2  gkg'1)  remained  almost  200m  above  the  boundary  layer  for  FS5.  There  is  also  greater 
moisture  in  the  boundary  layer  offshore  (9  gkg'1)  than  near  the  coast  (7.1  gkg'1) 
corresponding  to  higher  SST  away  from  the  coast.  At  the  top  of  the  boundary  layer,  there 
is  an  inversion  of  almost  10K  in  approximately  50m  depth.  Potential  temperature 
increases  to  approximately  305K  above  the  boundary  layer.  The  near-surface  0  are  about 

2K  colder  near  the  coastal  waters  than  offshore  (Figures  20,  21c)  and  consistent  with  the 
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water  temperature  below.  The  first  sounding  was  at  takeoff  from  the  Marina  airport  and 
experienced  slightly  higher  0  (about  2K)  then  the  upwelled  coastal  waters.  The  wind 
(Figure  2 Id)  is  strongest  offshore  (8  m  s'1)  in  the  boundary  layer,  decreases  to  2  m  s'1 
above  the  cloud  layer.  Each  of  the  variables  displayed  typical  stratocumulus  profiles  (e.g., 
Nicholls  1984). 
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Figure  21 .  Contour  plots  for  vertical  profiles  of  (a)  qc,  (b)  q,  (c)  6,  and  (d)  wind  speed 
for  RF02  full  soundings  (FS1-FS5)  and  short  soundings  (SSI)  combined. 
Flight  track  is  represented  by  solid  magenta  line. 

Black  line  indicates  slope  of  BL. 


2.  RF10  (26  September  2012) 

The  tenth  flight,  on  26  September  2012,  was  the  last  cloud  study  mission  in 
typical  summertime  solid  stratocumulus  cloud.  Figure  22  shows  three  consecutive 
GOES- 15  visible  satellite  images  of  Monterey  Bay.  Solid  cloud  was  throughout  the 
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duration  of  RF10.  Long  linear  transects  were  flown  both  in  cloud  and  near  the  surface  to 
get  a  full  boundary  layer  picture.  A  MODIS  image  at  2145  UTC,  with  RF10  flight  track 
overlaid  and  the  Full  soundings  color  coded,  is  presented  in  Figure  23.  Roll  structures  in 
cloud  aligned  with  the  predominant  northwesterly  wind  direction  are  apparent  in  all 
satellite  images.  West  of  the  flight  track  the  cloud  is  broken,  while  the  cloud  is  solid  for 
the  entire  flight  track  to  the  coast  and  clears  over  land. 


Figure  22.  Three  sequential  GOES-15  visible  images  on  26  September  2012,  at  1845, 
1945,  and  2030  UTC,  respectively,  show  cloud  variability  during  RF10.  White 
lines  are  the  flight  track.  Plane  icon  is  the  location  of  TO  at  time  of  image. 


Figure  23.  26  September  2012  MODIS  2145  UTC  image  overlaid  on  Google  Earth  with 

RF10  flight  track  (blue  solid  line)  and  FS1-FS13.  The  color  here  for  the  flight 
track  of  each  profde  is  the  same  as  shown  in  Figure  24. 
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Figure  24.  Same  as  Figure  19,  except  for  soundings  of  RF10. 


The  qc ,  9 ,  q ,  wind  speed  and  direction  for  each  of  the  13  full  soundings  are  plotted 
in  Figure  24.  The  vertical  profiles  of  qc  show  the  thickest  cloud  of  all  UPPEF  2012  cases 
with  maximum  cloud  thickness  approximately  330m.  The  highest  cloud  top  was  at  780m 
at  the  western  most  portion  of  flight  track  (leg  FS7),  and  the  lowest  cloud  top  was  at 
420m  near  the  coast  (FS12).  The  difference  in  cloud  top  heights  is  360m,  over  a  distance 
of  approximately  325km.  The  boundary  layer  height  plot  in  Figure  25,  though  based  on 
inversion  height,  represents  the  top  of  the  cloud  well,  since  it  was  solid  cloud  throughout 
the  flight  track.  The  descending  slope  of  the  boundary  layer  from  west  to  east  is  well 
characterized.  This  slope  is  quantified  using  contour  plots  discussed  below. 
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Figure  25.  Boundary  layer  heights  for  RF10.  Dotted  line  represents  flight  track. 

Coloreddots  and  numbers  represent  boundary  layer  height. 

All  of  the  full  soundings,  short  soundings,  and  level  legs  in  RF10  are  combined  to 
produce  the  contour  plots  in  Figure  26.  Figure  26a  shows  solid  cloud  from  all  soundings 
throughout  the  flight  and  the  solid  black  line  showing  the  slope  of  the  cloud  tops.  The 
same  solid  black  line  is  also  drawn  on  other  subplots  to  ensure  that  the  boundary  layer  top 
seen  from  all  variables  are  consistent.  Figure  26  shows  that  the  slope  out  to  -123.9W 
(200  km  offshore)  is  2  m  per  kilometer,  while  it  levels  off  further  west  at  about  0.5  m  per 
kilometer.  The  cloud  thickness  is  fairly  uniform  until  the  westernmost  part  of  the  track, 
beyond  which  point  the  cloud  becomes  broken.  The  full  soundings  did  not  extend  as  high 
as  in  RF02  to  allow  a  detailed  representation  of  the  free  atmosphere  above  the  BL,  but  the 
spatial  variability  is  well  shown  in  the  BL.  Moisture  variability  across  the  boundary  layer 
top  varies  substantially  as  seen  in  Figure  26b.  The  BL  is  moist  at  9  gkg'1,  but  the  free 
atmosphere  above  is  not  as  dry  as  in  RF02,  resulting  in  much  smaller  moisture  jump 
across  the  cloud  top.  At  the  coast,  q  decreases  more  than  over  the  ocean.  Moisture  is 
mixed  above  the  BL  and  for  most  soundings,  q  is  greater  than  5  gkg'1.  In  comparison,  the 
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q  values  were  all  less  than  5  gkg'1  above  the  BL  during  RF02.  A  vertical  profile  of  0  is 
shown  in  Figure  26c.  The  slope  of  the  potential  temperature  inversion  is  consistent  with 
the  slope  of  the  cloud  top.  In  the  BL,  0  is  uniform  with  height,  with  a  strong  inversion  of 
about  10K  at  the  BL  top.  Above  the  BL,  the  0  increases  to  300-305K.  The  darker  blue  is 
lower  values  of  0  and  the  red  colors  are  higher  values  of  0.  Finally,  Figure  26d  is  of  wind 
speed.  As  in  RF02,  the  wind  speeds  are  typically  higher  in  the  BL  and  decreasing  above 
by  2-3  ms'1.  The  wind  is  weakest  near  the  coast,  at  approximately  4-6  ms'1.  Offshore,  the 
winds  increase  to  10-12  ms'1.  The  slope  of  wind  shear  is  fairly  consistent  with  the  top  of 
the  cloud  layer,  except  for  at  the  westernmost  portion  of  the  track,  where  the  wind 
increases  slightly  above  the  cloud  top. 
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Figure  26.  Contour  plots  for  vertical  profiles  of  (a)  qc.  (b)  q,  (c)  0,  and  (d)  wind 
speed  for  entire  RF10  flight  track.  Flight  track  is  represented  by 
solid  magenta  line. 


47 


B.  SST  VARIABILITY  AND  ITS  IMPACT  ON  LOW  LEVEL  TURBULENCE 

During  each  of  the  flights,  SST  data  was  collected  from  the  downward  looking 
infrared  radiometer.  The  SST  data  from  the  lowest  flight  levels  (below  100m)  of  the  first 
10  flights  are  shown  in  Figure  27.  The  range  of  the  SST  in  all  figures  is  the  same  as 
shown  in  the  color  bar  next  to  SST  from  RF10. 
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Figure  27.  RF01-RF10  contour  plot  for  SST.  Solid  magenta  lines  represent  level  legs 

at  an  altitude  of  100m  or  less.  SST  is  in  °C. 
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Figure  27  shows  significant  spatial  and  temporal  SST  variability  in  the  UPPEF 
2012  measurement  area.  The  frequent  presence  of  SST  front  is  quite  apparent  as  seen  in 
nearly  all  flights  except  for  RF02.  The  magnitude  of  variability  in  SST  is  in  general  1- 
2.5°C.  During  the  first  flight,  RF01,  there  was  a  strong  SST  front  near  (35.85N 
123.45W).  The  coldest  and  rather  uniform  SST  field  was  observed  on  RF02  under  solid 
stratocumulus  clouds.  There  was  a  weak  SST  front  of  a  1°C  change,  near  (36. 3N 
122. 9W),  during  RF03.  Another  weak  front  was  evident  in  RF04,  near  (36. IN  122. 7W). 
During  RF05,  the  northern  leg  was  colder  than  the  southern  leg  in  the  U/box-pattem  by 
approximately  2°C.  The  sixth  flight  exhibited  a  2°C  SST  change  along  the  surface  tracks, 
but  no  sharp  front  was  observed.  RF07  was  primarily  in  relatively  warm  SST,  near  16°C. 
The  middle  of  RF08  had  SST  greater  than  16°C,  surrounded  by  1°C  colder  SST  to  the 
southwest  and  2°C  colder  SST  to  the  northeast.  There  is  a  weak  SST  front  with  a  1°C 
change  in  SST  during  RF09  along  the  eastern  part  of  the  U-pattem,  near  122. 9W.  During 
RF10,  a  strong  SST  front  with  a  2.5°C  change  in  SST  near  (36.55N  123. 7W)  was 
observed.  The  sharp  SST  front  during  RF01  and  its  impact  on  near-surface  mean  and 
turbulence  are  examined  in  detail,  followed  with  similar  discussions  on  SST  fronts  during 
RF03  and  RF05. 
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1. 


Strong  SST  Front  in  RF01 


Aug  31  2012  US  West  Coast  JPL  G1SST 


Figure  28.  NASA  JPL  1km  Global  SST  (G1SST)  for  31  August  2012  overlaid  on 
Google  earth  with  RF01  flight  track  (light  blue  line),  with  LL7  and  LL12 
(purple  lines),  LL8  and  LL9  (orange  lines),  and  LL10  and  LL1 1  (red  lines). 


A  graphical  representation  of  the  Global  1  km  SST  (G1SST)  field  provided  by  the 
National  Aeronautics  and  Space  Administration  (NASA)  Jet  Propulsion  Laboratory  (JPL) 
at  California  Institute  of  Technology’s  Regional  Ocean  Modeling  System  (ROMS)  group 
for  the  gap-free  blended  SST  along  the  US  West  Coast  was  downloaded  for  each  UPPEF 
2012  flight.  The  image  for  31  August  2012  is  overlaid  in  Google  Earth  with  the  RF01 
flight  track  and  level  legs  in  the  cold  and  warm  sectors  of  the  SST  front  (Figure  28).  The 
sharp  SST  front  is  evident  in  the  image  along  the  red  and  purple  legs  and  corresponds  to 
the  SST  front  revealed  by  the  UPPEF  2012  data  (Figure  29).  Two  levels  of  legs  were 
flown  over  the  SST  front,  LL9  and  LL10  were  flown  just  under  100m,  while  LL7,  LL8, 
LL1 1,  and  LL12  were  flown  closest  to  the  surface  (~30m).  Time  series  of  SST,  0,  q,  wind 
speed,  wind  direction,  and  pressure  comparing  the  four  legs  are  shown  in  Figure  30. 
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Figure  29.  Contour  plot  for  SST  during  RF01.  Flight  tracks  are  represented  by  magenta 
solid  line.  Near  surface  LL7,  LL12,  LL10  and  LL1 1  crossed  over  the  strong 
SST  gradient.  Near  surface  LL8  and  LL9  are  in  the  warm  SST  region,  as 

indicated. 
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Figure  30.  Time  series  of  SST,  6 ,  q ,  wind  speed,  wind  direction,  and  pressure  plotted 
against  longitude  for  (a)  LL7  and  LL12,  and  (b)  LL10  and  LL1 1  crossing 

the  SST  front. 
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Figure  3 1 .  Same  as  in  Figure  30,  except  for  LL8  and  LL9  in  the  warm  SST  region. 

The  time  series  plots  are  used  to  compare  legs  to  evaluate  temporal  and  spatial 
variability  of  the  SST  front  during  RF01  for  LL7  and  LL12  (Figure  30a),  LL10  and  LL1 1 
(Figure  30b),  and  LL8  and  LL9  (Figure  31).  LL7  and  LL12  overlay  each  other  at  an 
altitude  of  35m,  with  LL7  passing  the  SST  front  53  minutes  before  LL12.  LL10  was  at 
93m  and  20  minutes  before  LL1 1  at  38m  in  altitude.  Overall,  the  time  series  match  fairly 
well  for  LL7  and  LL12,  and  suggests  not  much  temporal  variability  for  the  SST  front 
during  the  flight,  except  the  wind  direction  had  a  more  northerly  component  in  LL7.  For 
LL10  and  LL11,  all  variables  from  the  warm  sector  showed  similar  magnitude  at 
different  height,  suggesting  that  the  boundary  layer  is  well  mixed  to  ~100m  over  the 
warm  water.  This  is  consistent  with  the  results  shown  in  Figure  31,  where  the  entire  legs 
at  both  levels  were  from  the  warm  region.  Figure  30b  indicates  apparent  stratifications  in 
potential  temperature  (0),  water  vapor  (q),  and  wind  speed  over  the  cold  water.  The 
section  over  cold  water  from  the  low  level  (LL1 1)  shows  larger  specific  humidity,  colder 
temperature,  weaker  winds  compared  to  the  respective  variables  from  the  layer  above,  all 
connected  with  the  cooler  surface. 
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The  response  of  boundary  layer  turbulence  to  the  significant  SST  variability  is 
discussed  here.  In  all  of  the  plots,  turbulence  is  greater  in  the  warm  SST  region  (west  of 
front)  compared  to  the  cold  SST  regions  (east  of  front).  The  sensible  heat  flux  (SHF), 
latent  heat  flux  (LHF),  wind  stress,  difference  of  T  and  SST,  and  SST  for  LL7  and  LL12 
are  shown  in  Figure  32.  Similar  plots  were  generated  for  LL10  and  11  in  Figure  33  and 
for  LL8  and  9  in  Figure  34.  Table  3  provides  the  calculated  mean  values  of  each  variable 
for  all  relevant  legs  prepared  into  warm  (top  table)  and  cold  (bottom  table)  regions.  A 
quick  glance  at  Table  3  indicates  significant  differences  in  all  variables  similar  to  those 
shown  in  Figures  30-34.  The  warm  region  has  near  surface  (~30  m)  SHF  in  the  range  of 
6-13  Wm'2,  LHF  of  90-118  Wm'2,  and  stress  between  0.30-0.37  Nm'2,  and  negative  air- 
SST  difference  (unstable  stratification).  In  contrast,  the  cold  water  region  are  stably 
stratified  with  air-sea  temperature  difference  between  0.6  and  1.8°C,  resulting  in  small 
but  consistently  negative  SHF  at  both  levels.  The  wind  stress  is  very  weak  over  the  cold 
water  region,  too.  Latent  heat  flux  are  small,  too,  but  consistently  positive.  Water  vapor 
supply  to  the  boundary  layer  over  the  cold  region  is  definitely  much  weaker  compared  to 
the  warm  region. 

Wind  speed  was  observed  to  be  consistently  weaker  over  the  cold  region 
compared  to  its  warm  neighbor  (Table  3).  The  weaker  speed  on  the  cold  side  of  a  SST 
front  is  primarily  due  to  the  stable  boundary  layer  formed  over  cold  SST,  which  leads  to  a 
very  shallow  BL  and  thus  a  strong  momentum  flux  convergence  to  reduce  the  wind  speed 
(Skyllingstad  et  al.,  2007).  Some  flux  divergence  are  seen  over  the  warm  water  region, 
the  divergence  is  not  apparent  over  the  cold  water  region,  likely  imbedded  in  the 
sampling  error  in  the  fluxes  of  smaller  magnitudes.  . 
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Figure  32.  Spatial  variability  showing  SHF,  LHF,  wind  stress,  difference  of  air 

temperature  and  SST  (T-SST),  and  SST  for  (a)  LL7  and  (b)  LL12  crossing 

the  SST  front  during  RF01. 


Figure  33  shows  the  flux  variations  of  air  temperature  from  SST,  and  SST  based 
on  longitude  of  the  overlying  LL10  and  11.  The  main  difference  between  this  pair  and 
LL7  and  12  (Figure  32)  is  the  height  difference  in  the  level  legs.  LL7  and  12  were  at  the 
same  altitude,  while  LL10  is  55m  above  LL11.  The  SST  front  is  located  at  123.4W  in 
Figure  33,  consistent  with  Figure  32.  The  SHF  is  almost  double  right  before  the  SST 
front  in  the  warm  sector  closer  to  the  surface  in  LL1 1,  at  123. 5W,  than  at  93m  for  LL10. 
The  mean  values  show  that  the  SHF  values  decrease  with  altitude,  while  the  LHF 
remain  approximately  the  same  (Table  3).  The  u-wind  stress  is  also  more  negative  from 
123.4-123. 5W  in  LL1 1,  than  LL10  as  it  is  approximately  zero.  The  u  and  v  wind  stress  at 
the  surface  nearly  reflect  each  other  about  the  axis  for  LL11,  while  the  wind  stress  for 
LL10  is  not  symmetric.  This  wind  stress  difference  is  not  as  evident  in  the  time  series 
plots  (Figure  30b)  as  in  the  flux  plots.  The  wind  backs  by  3°  as  both  legs  pass  the  SST 
front  form  warm  to  cold,  and  decreases  in  speed  by  1.3°C  in  LL9  and  only  0.3°C  in  LL1 1 
(Table  3). 
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Figure  33.  Same  as  in  Figure  32,  except  for  (a)  LL10  and  (b)  LL1 1  crossing  the 

SST  front  during  RF01. 


Figure  34  displays  the  same  plots  as  Figures  32  and  33,  but  for  LL8  and  9.  LL8 
and  9  were  entirely  located  in  the  warm  SST  region.  LL9  was  at  90m  in  altitude,  which 
was  59m  above  LL8,  at  31m  in  altitude.  Both  show  positive  SHF  and  LHF  fluxes,  but 
consistent  with  LL10  and  11,  both  fluxes  decrease  with  altitude.  The  LHF  for  LL9  is 
about  half  of  that  for  LL8  and  the  SHF  is  about  a  third  less  for  LL9  than  LL8  (for  exact 
values,  see  Table  3).  The  wind  stress  does  not  change  much,  as  the  wind  direction  and 
speed  are  consistent  for  both  legs.  The  temperature  difference  is  almost  0.5°C  and 
0.5K  more  at  90m  than  31m,  consistent  with  LL10  and  LL11.  The  wind  speed  was 
stronger  for  LL8  by  at  least  0.5  ms"1  than  all  other  legs  (Table  3).  The  wind  direction  is 
consistent  in  the  warm  SST  for  LL8,  LL9,  LL7,  and  LL12. 
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Figure  34.  Same  as  in  Figure  32,  except  for  (a)  LL8  and  (b)  LL9  in  the  warm  SST  region 

during  RF01. 


Table  3.  Comparison  of  mean  and  turbulent  fluxes  over  warm  and  cold  SST  regions. 


RF01  Mean 

variables 

LL11  (z=38m) 
Warm 

LL10  (z=93m) 
Warm 

LL7  (z=35m) 
Warm 

LL12  (z=35m) 
Warm 

LL8  (z=31m) 
Warm 

LL9  (z=90m) 
Warm 

SHF  (Wm’2| 

9.0964 

6.4384 

10.2287 

13.2936 

5.8609 

4.1867 

LHF  (WrrT2) 

117.8148 

117.1868 

104.0729 

89.4846 

100.1591 

43.9612 

Stress  (NrrT2) 

0.3621 

0.3592 

0.2957 

0.3741 

0.1022 

0.1061 

Wind  Speed  (ms’1) 

8.933 

9.3587 

8.9654 

8.5702 

9.9017 

9.388 

Wind  Direction  (°) 

327.4332 

320.9001 

329.9618 

320.5 

329.6294 

330.5528 

e-SST(K) 

-0.4932 

-1.0551 

-0.6582 

-0.471 

-0.3841 

-0.8582 

SST  (°C) 

15.9801 

15.9686 

16.0799 

15.9975 

15.9687 

15.8803 

RF01  Mean 

variables 

LLll  (z=38m) 
Cold 

LL10  (z=93m) 
Cold 

LL7  (z=35m) 
Cold 

LL12  (z=35m) 
Cold 

SHF  (WrrT2) 

-0.9813 

-0.6611 

-1.5665 

-1.3253 

LHF  (WrrT2) 

3.263 

4.5991 

4.7319 

5.227 

Stress  (Nm'2) 

0.0268 

0.0319 

0.024 

0.0352 

Wind  Speed  (ms’1) 

7.4146 

9.0964 

7.6675 

7.2658 

Wind  Direction  (°) 

320.786 

317.434 

325.0623 

312.0479 

0  -SST  (K) 

0.9381 

0.6315 

1.6552 

1.8199 

SST  (°C) 

13.7164 

13.758 

13.0769 

12.977 
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2. 


RF03  -  Weaker  and  Warmer  SST  Front 


The  third  flight,  RF03,  is  another  example  of  measurements  over  a  weaker  SST 
front.  The  SST  remains  above  14°C  in  the  cold  region,  warmer  than  that  found  in  RF01. 
The  G1SST  image  overlaid  on  Google  earth  with  RF03  flight  track  and  front  passing  legs 
(LL16  and  LL20)  highlighted  is  in  Figure  35.  The  SST  variability  shown  in  the  G1SST 
data  does  not  match  well  with  the  data  collected  by  the  TO  in  Figure  36.  It  appears  to  be 
about  1-2°C  too  low  and  the  warm  eddy  of  SST  is  too  far  to  the  west  for  the  G1SST.  This 
is  likely  associated  with  the  errors  in  the  satellite  data.  The  cloud  cover  for  the  third  flight 
was  much  more  extensive  than  for  RF01,  which  plays  a  factor  in  the  accuracy  of  satellite 
derived  SST.  Only  two  level  legs  collected  data  from  this  SST  front,  LL16  at  35m  and 
LL20  at  92m  (Figure  36).  There  was  a  2°C  change  at  SST  front,  from  15°C  to  17°C 
(Figure  36).  Spatial  variations  of  SST,  0,  q,  wind  speed,  and  wind  direction  based  on  the 
longitude  of  LL16  and  LL20  are  also  shown  in  Figure  37.  Corresponding  flux  and  flux- 
related  quantities  from  LL16  and  LL20  are  presented  in  Figure  38. 


Figure  35.  NASA  JPL  G1SST  for  06  September  2012  overlaid  on  Google  earth  with 
RF03  flight  track  (light  blue  line),  LL16  (green  line),  and  LL20  (red  line). 
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RF03  COMB-SST 
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Figure  36.  Contour  SST  plot  for  RF03.  Flight  track  is  the  magenta  solid  line.  LL16  and 
LL20  overlaid  each  other  and  crossed  over  the  strongest  SST  gradient  as 

indicated. 
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Figure  37.  Same  as  in  Figure  30,  except  for  LL16  and  LL20  of  RF03. 
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Figure  38.  Same  as  in  Figure  30,  except  for  (a)  LL16  and  (b)  LL20  crossing  the  SST 

front  during  RF03. 

Figure  36  shows  the  location  of  the  SST  front  at  (36. 3N  122. 9W)  during  RF03. 
Compared  to  RF01,  the  results  from  RF03  differ  significantly.  The  values  for  SHF  and 
LHF  are  positive  for  the  entire  leg  for  LL16  at  an  altitude  of  36m,  vice  just  in  the  warmer 
SST  region,  shown  in  the  results  for  RF01.  The  SHF  is  slightly  negative  in  the  warm  SST 
region  for  LL20  at  92m  in  height.  Once  the  SST  front  is  passed  in  the  colder  SST,  the 
SHF  becomes  positive — opposite  the  results  of  RF01.  This  is  likely  due  to  LL20  being 
above  the  surface  layer.  The  surface  layer  is  generally  10%  of  the  BL  height.  The  BL 
height  is  approximately  500m  over  LL20,  making  the  surface  layer  at  50m  in  height.  The 
level  leg  height  of  90m  is  well  above  the  surface  layer;  other  factors  than  surface  forcing 
is  contributing  to  the  negative  flux.  Also,  another  difference  from  RF01  is  that  the  air 
temperature  difference  from  the  SST  is  negative  on  both  sides  of  the  front.  The  SST  is 
warmer  than  the  overlying  air  for  both  legs,  thus  should  always  producing  some  upward 
vertical  motion  and  heat  fluxes.  A  third  difference  from  RF01  is  that  the  wind  speed  is 
rather  light  for  RF03,  less  than  5  ms'1  during  both  LL16  and  20,  while  the  wind  speed 
during  RF01  was  almost  double.  This  light  wind  speed  shows  in  the  nearly  zero  wind 
stress,  contributing  to  the  low  fluxes. 
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3. 


RF05  -  Gradual  SST  Change  and  Cold  Eddy 


Another  interesting  flight  showing  SST  variability  was  RF05.  Clouds  were 
primarily  solid  to  broken  on  11  September  2012.  The  SST  variability  shown  in  the 
G1SST  image  matches  fairly  well  with  the  TO  data  (Figures  39  and  40).  There  was  not  a 
single  distinct  SST  front  during  this  flight,  but  a  gradual  increase  in  SST  to  15.3°C  from 
the  east  in  the  southern  legs  and  the  northern  legs  of  the  box  appear  to  cross  a  filament  of 
13°C  (Figure  40).  There  were  four  legs  that  collected  data  near  the  surface  showing  SST 
variability  during  RF05.  LL9  was  at  89m  and  overtop  LL10  at  35m  in  altitude.  LL1 1  and 
LL12  were  at  33m  and  35m,  respectively,  and  were  two  other  surface  legs  making  up  the 
box  pattern  flown,  as  shown  in  Figure  40. 


Sep  11  2012  US  West  Coast  JPL  G1SST 


Figure  39.  NASA  JPL  1km  Global  SST  for  1 1  September  2012  overlaid  on  Google  earth 
with  RF05  flight  track  (light  blue  line),  with  LL10,  LL1 1  and  LL12 
(green  lines),  LL9  and  LL13  (red  lines). 


60 


RF05  COMB-SST 


Figure  40.  Contour  plot  of  RF05  SST  with  LL9,  LL10,  LL1 1,  LL12  and 

height  indicated  on  the  track. 


Figure  41  shows  the  time  series  of  SST,  0,  q,  wind  speed,  wind  direction,  and 
pressure  for  the  overlying  pair  (LL9  and  LL10),  LL1 1,  and  LL12.  LL9  and  LL10  show  a 
gradual  decrease  in  SST  to  the  east  towards  the  coast.  LL11  shows  a  gradual  decrease  in 
SST  to  the  northwest,  possibly  entering  a  cold  filament.  Wind  is  northwesterly  between 
5-10  ms'1  for  all  four  legs,  with  the  wind  speed  at  LL9  slightly  higher  than  the  other  legs. 
Wind  direction  is  consistent  for  all  four  legs.  SST  drops  by  about  1°C  in  a  small  portion 
of  LL12  from  122.75W-122.95W  from  the  surrounding  SST.  Turbulence  is  evident  in 
SST  greater  than  approximately  15°C.  An  interesting  feature  for  LL9  and  LL10,  is  that  as 
the  SST  decreases,  0,  q,  and  wind  speed  also  decrease.  These  variables  contribute  to  the 
SHF,  LHF,  and  wind  stress  profdes  for  each  leg  in  Figure  42. 
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Figure  41 .  Same  as  in  Figure  30,  except  for  (a)  LL9-LL10,  (b)  LL1 1,  and 

(c)  LL12  of  RF05. 


62 


RFOB  201  20911  LL9  RFQ5  201  2091 1  LL ID  2=35 


-122.9  -122.8  -122,7  -122.9  -122.5  -122,4  -122.3  -122.2  -122.1  -123  -122,8  -122.9  -122,4  -122.2  -122 

Longitude  (°)  Longitude 


RFQ5  201  20911  LL1 1  i-33  RF05  29120911  LL12  3-36 


-123.1  -123.05  -123  -122.95  -122.9  -122.85  -123.2  -123  -122.8  -122.6  -122.4  -122,2 

Long  Itude  {“*)  Long  itude  {*) 


Figure  42.  SHF,  LHF,  wind  stress,  temperature  and  SST  difference,  and 

SST  profiles  for  (a)  LL9,  (b)  LL10,  (c)  LL11,  and  (d)  LL12  during  RF05. 


The  SHF  variation  for  all  four  legs  is  consistent  with  the  results  from  RF01.  When 
the  temperature  difference  is  negative  (unstable  stratification),  there  are  positive  SHF  and 
LHF.  When  the  difference  is  positive,  then  there  are  no  fluxes  or  they  are  slightly 
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negative.  The  LHF  values  are  very  low  compared  to  RF01,  along  with  the  wind  stress 
being  at  a  minimum.  There  are  not  many  differences  between  LL9  and  LL10,  except  that 
there  is  slightly  greater  wind  stress  in  LL9  than  LL10.  The  results  from  LL9  at  89m  has 
the  same  trend  as  the  other  legs  and  RF01,  indicating  that  this  leg  is  possibly  in  the 
surface  layer.  The  BL  height  near  LL9  was  at  800m.  Using  the  general  rule  of  thumb  for 
the  surface  layer  height  at  10%  of  the  BLH,  the  results  suggest  that  the  surface  layer  may 
extend  up  to  90m  in  height. 

C.  TURBULENCE  CHARACTERISTICS  IN  CLOUD  AND  SURFACE 

The  second  flight  is  a  prime  example  with  sufficient  level  legs  to  get  an  accurate 
depiction  of  turbulent  statistics  in  solid  cloud  and  its  associated  boundary  layer.  The 
mean  data  in  each  level  leg  is  represented  by  a  circle  in  the  flux  and  variance  profiles  in 
Figures  43  and  44.  The  vertical  axes  are  mean  leg  height  scaled  by  BLH  at  the  center 
location  of  the  leg  .  Figure  43a  is  the  sensible  heat  flux  profile,  w'9' ,  and  proportional  to 
the  buoyancy  flux.  Figure  43b  is  the  moisture  flux  profile,  wrqr,  Figure  43a  shows  that 
the  largest  buoyancy  flux  is  in  the  solid  cloud,  indicating  that  radiative  cooling  in  the 
cloud  top  is  the  primary  driver  of  turbulence  in  the  STBL.  The  surface  layer  has  a  small 
positive  buoyancy  flux  over  cold  water.  Compared  to  the  solid  cloud,  the  surface  is  not  a 
significant  contributor  to  turbulence  in  the  STBL.  However,  the  surface  is  the  primary 
source  of  moisture  into  the  STBL,  as  represented  by  the  large  surface  moisture  flux  in 
Figure  43b.  This  figure  also  shows  that  there  is  a  significant  amount  of  water  vapor  flux 
in  solid  cloud.  The  variance  profiles  for  w  and  9  are  shown  in  Figure  44a  and  44b, 
respectively.  Strong  turbulent  eddies  in  the  cloud  layer  are  shown  in  Figure  44a.  The 
larger  temperature  perturbations  in  Figure  44b,  are  in  the  cloud  layer.  The  cold 
downdrafts  originating  at  the  cloud  top  contribute  to  large  variance.  A  scatter  plot  is  used 
to  examine  how  the  individual  perturbations  contribute  to  the  turbulence  characteristics. 
Figure  45  represents  the  relationship  between  w’  and  9’  in  solid  cloud  from  an  in-cloud 
leg  (LL1)  during  RF02.  The  results  are  oriented  in  the  first  and  third  quadrants,  so  there 
are  either  warm  updrafts  or  cold  downdrafts.  The  scatter  plot,  along  with  the  individual 
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probability  distribution  functions  (PDF),  however,  is  strongly  negatively  skewed  and 
indicates  that  the  cold  downdrafts  are  the  primary  forcing  of  turbulence  in  the  cloud. 
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Figure  43.  Flux  profdes  of  (a)  w'9'  and  (b)  w'q'  for  RF02.  The  vertical  axis  is  height 

scaled  by  BLH. 


Variance  Profiles  RF02 


Figure  44.  Same  as  in  Figure  43,  except  for  variance  profiles  of  (a)  w'2  and  (b)  9'2 

for  RF02. 

The  difference  between  the  contribution  from  warm  SST  and  solid  cloud  is 
significant  and  is  shown  in  Figure  46.  The  joint  PDF  for  RF02  in  solid  cloud  is  compared 
to  RF01  over  the  warm  SST.  Figure  46a  shows  the  PDF  for  vv  ’  and  9  ’  and  Figure  46c  is 
for  w  ’  and  q  ’  in  the  warm  SST  region  of  LL1 1  (LL1 1 W)  in  RF01 .  Figure  46b  shows  the 
PDF  for  w’  and  9’  and  Figure  46d  for  w’  and  q ’  in  the  solid  stratocumulus  cloud  during 

LL1  in  RF02.  Comparing  the  PDFs,  it  can  be  easily  seen  that  those  representing  the  solid 
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cloud  are  negatively  skewed  (Figure  46b,  d)  and  the  PDFs  for  the  legs  above  the 
relatively  warm  SST  are  positively  skewed  (Figure  46a,  c).  The  positive  skewness 
indicated  warm  and  moist  updrafts  from  the  surface  over  relatively  warm  SST.  The 
negative  skewness  indicated  cool  and  dry  downdrafts  in  the  cloud. 


RF02  20120904  LL1  z=447m 


w’ 


Figure  45.  Scatter  plot  for  solid  cloud  leg  (LL1)  of  w  ’  and  6  ’  for  RF02, 

showing  negative  skewness. 
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Figure  46.  Joint  PDF  for  comparing  (a)  w’  and  9’  for  LL1 1W  in  RF01,  (b)  w’  and  9’  for 
LL1  in  RF02,  (c)  w  ’  and  q  ’  for  LL1 1 W  in  RF01,  and  (d)  w  ’  and  q  ’  for  LL1  in 

RF02. 
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V.  DISCUSSION 


A.  SUMMARY  AND  CONCLUSIONS 

A  major  objective  of  this  thesis  was  to  provide  the  initial  analysis  of  the  UPPEF 
2012  aircraft  measurements  by  examining  all  sampled  cases  and  identify  specific  cases 
for  future  in-depth  study.  In  this  thesis,  the  boundary  layer  heights  were  estimated  for  all 
flights  which  show  a  descending  slope  towards  the  coast.  Two  of  the  flights,  RF02  and 
RF10,  were  analyzed  in  depth.  Results  from  these  two  flights  show  that  the  BL  decreases 
at  a  rate  of  about  2m  km'1  to  the  coast  from  approximately  200km  offshore.  Westward  of 
200km,  the  slope  levels  out  to  a  rate  of  approximately  0.5  m  km'1,  as  seen  in  the  case 
study  of  RF10.  Conceptually  and  qualitatively,  the  boundary  layer  is  easily  understood, 
but  accurately  quantifying  the  numerous  physical  processes  that  affect  the  BL  is  much 
harder  to  do.  Next,  the  sea  surface  temperature  variability  for  each  flight  was  plotted  and 
three  flights  were  chosen  to  analyze  in  detail.  There  was  a  strong  SST  front  during  RF01, 
a  weaker  and  warmer  front  during  RF03,  with  a  gradual  change  in  SST  during  RF05. 
Each  of  the  surface  legs  covering  the  changes  in  SST  were  evaluated  and  compared  using 
the  spatial  variability  of  relevant  variables  at  different  levels.  This  research  found  that  the 
small  scale  variability  in  the  SST  resulted  in  significant  variation  in  the  surface  exchange 
of  sensible  and  latent  heat  fluxes,  and  wind  stress.  Such  variability  makes  it  difficult  to 
correctly  parameterize  surface  fluxes.  The  presence  of  cool  downdrafts  in  the  upper  cloud 
layer  is  evident  in  the  joint  probability  density  distribution  and  buoyancy  flux  profiles. 
The  results  from  this  research  will  guide  future  cloud  parameterization  scheme 
development  in  forecast  models. 

B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

The  findings  of  this  initial  study  yield  intriguing  results  of  the  boundary  layer  and 
its  characteristics  under  solid  to  broken  stratocumulus  clouds  and  in  clear  conditions. 
There  is  a  lot  of  work  to  be  done  to  further  understand  the  extensive  amount  of  data 
collected  during  UPPEF  2012.  One  such  area  is  in  quantifying  the  entrainment  events 
frequently  seen  during  the  initial  analyses.  There  are  numerous  examples  in  the  data  that 
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show  breaks  in  the  LWC,  corresponding  to  breaks  in  the  number  of  concentration,  and 
decrease  in  q,  amongst  others.  These  breaks  are  possibly  due  to  the  dry  downdrafts 
resulting  from  entrainment.  The  EIL  in  the  solid  cloud  needs  to  be  identified,  along  with 
cloud  top  jumps.  Also,  there  are  many  interesting  features  in  the  cloud  from  the  satellite 
images.  These  structures  need  to  be  analyzed  further  to  reveal  scales  of  variability  in  the 
cloud.  Finally,  the  surface  flux  parameterization  over  an  inhomogeneous  ocean  surface 
needs  to  be  examined  in  further  detail. 
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